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SUMMARY 


The  general  problem  of  aircraft  electrostatic 
charging  and  one  of  the  more  efficient  discharger  systems 
presently  available  are  discussed.  Methods  of  discharger 
system  design  and  analysis  are  developed  and  stated  for 
this  system.  These  methods  are  incorporated  into  design 
curves  and  other  data  to  facilitate  and  optimize  design*. 

The  prcolem  of  efficient  discharge  corona  electrode 
design  is  treated  extensively.  With  the  information  on 
corona  developed  in  the  report,  it  becomes  possible  to 
include  the  corona  electrode  in  the  system  design  using 
conventional  engineering  design  methods >  Data  permitting 
optimization  of  corona  electrode  performance,  that  is, 
current  maximization  at  a  given  electrode  voltage  is  given. 
The  effects  of  aircraft  altitude  are  significant  on  electrode 
operation. 

Numerous  approaches  to  the  reduction  of  system 
voltage  requirements  are  discussed. 
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FOREWORD 


The  Dyriasciences  Corporation,  under  the  sponsorship 
of  the  U.S.  Army  Aviation  Materiel  Laboratories*  (USAAML), 
Fort  Eustis,  Virginia,  has  conducted  a  study  to  provide  the 
theory  for  the  efficient  design  and  optimization  of  aircraft 
electrostatic  discharging  systems.  The  investigators  were 
Messrs.  Juan  de  la  Cierva,  Herbert  J.  Gillis,  and  Paul  B. 
Wilson,  Jr.  The  cognizant  USAAVLABS  project  officer  was  Mr. 
S.  Blair  Poteate,  Jr. 

The  project  was  conducted  during  the  period  June 
1963  through  January  1964. 

The  continuing  interest  of  the  U.S.  Army  Aviation 
Materiel  Laboratories  in  helicopter  performance  impro\ ement 
made  possible  this  contribution  to  aircraft  electrostatic 
discharging  and  is  gratefully  acknowledged. 


^Formerly,  U.S.  Army  Transportation  Research  Command. 
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SECTION  A:  INTRODUCTION 


As  the  role  of  the  helicopter  as  an  operational 
vehicle  expands,  it  becomes  increasingly  important  to  deal 
effectively  with  the  problem  of  large  electrostatic  charge 
accumulation  on  the  aircraft.  This  charge  is  primarily 
caused  by  triboelectric  effects  between  the  rotor  blades 
and  the  air  through  which  they  move.  Triboelectricity 
has  been  defined  as  the  electric  charge  transfer  which 
takes  place  between  dissimilar  substances  when  brought 
into  contact  and  subsequently  separated.  The  charge 
transfer  is  due  to  the  different  number  of  atomic  outer 
shell  electrons  which  migrate  from  one  substance  to  the 
other  and  vice  versa.  The  interception  of  solid,  liquid, 
or  gaseous  particles  in  the  air  by  the  helicopter  rotor 
blades  and  fuselage  produces  a  cumulative  triboelectric 
interchange  which  results  in  an  electrostatic  charge 
potential  and  electrical  energy  buildup  on  the  aircraft. 
The  large  energy  transfer  which  can  occur  when  the  charged 
aircraft  comes  into  contact  with  ground  may  present  a 
serious  hazard.  The  danger  of  electrical  shock  to  per¬ 
sonnel,  ignition  of  flammable  gas  mixtures,  and  detonation 
of  explosives  severely  affects  cargo  handling  and  per¬ 
sonnel  rescue  operations.  The  hazards  of  electrostatic 
charge  are  eliminated  if  the  stored  electrical  energy  of 
the  aircraft  is  kept  at  a  level  of  one  mill! joule  or  less 
(Reference  21,  page  7). 

Many  attempts  to  build  aircraft  electrostatic 
discharging  systems  have  been  made.  The  obvious  expedient 
of  discharging  the  aircraft  by  electrically  contacting  the 
earth  or  sea  by  cable  has  an  associated  spark  hazard  and 
introduces  handling  and  operating  problems.  It  is  there¬ 
fore  an  undesirable  solution  in  many  cases.  Another 
method  of  preventing  charge  accumulation  1 s  the  applica¬ 
tion  of  coatings  to  match  the  dielectric  constant  of  the 
aircraft  to  that  of  the  atmosphere  (Reference  2).  This 
method,  however,  has  proven  inadequate  in  preventing 
triboelectric  charging  in  all  operating  conditions  due  to 
the  fact  that  triboelectric  charging  is  intrinsically  a 
surface  pheonomenon,  and  triboelectric  surface  character- 
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istics  of  operational  aircraft  are  highly  dependent  on 
surface  conditions  such  as  dust,  oil,  grit  and  other 
contaminant  deposits,  as  well  as  snow  and  ice  accumula¬ 
tions.  These  conditions  vary  in  very  significant  amounts 
during  the  normal  use,  maintenance,  and  service  of  mili¬ 
tary  and  commercial  aircraft.  The  use  of  radioactive 
coatings  to  discharge  aircraft  by  rendering  the  adjacent 
atmosphere  conductive  because  of  ionization  has  also  been 
investigated  (Reference  2).  It  was  found  that  an  ex¬ 
cessively  large  volume  of  atmosphere  (50,000  cubic  meters) 
had  to  be  irradiated  if  the  maximum  health  safety  limit 
was  not  to  be  exceeded. 

The  use  of  passive  corona  discharger  devices, 
such  as  the  AN/ASA- 3  wick  dischargers  commonly  employed 
in  fixed  wing  aircraft  and  similar  devices,  has  been 
investigated  and  proven  to  be  impractical  for  helicopters 
(Reference  3).  The  high  level  of  aircraft  potential  re¬ 
quired  to  discharge  the  relatively  large  charging  currents 
frequently  encountered  in  helicopter  operations  by  means 
of  these  devices,  results  in  high  electrostatic  energies 
stored  on  the  aircraft  which  may  exceed,  by  several  orders 
of  magnitude,  the  safe  level  required. 

A  fully  successful  discharger  system  was  built 
by  the  Dynasciences  Corporation  under  USAAML  Contract 
DA-177-AMC-114(T) .  This  system  automatically  controls 
the  voltage  difference  between  a  corona  electrode  and  the 
aircraft  surface  as  a  function  of  the  electrostatic  field 
existing  on  the  aircraft’s  outer  surface.  The  electro¬ 
static  field  is  measured  by  a  sensitive  device  capable  of 
detecting  field  intensities  corresponding  to  aircraft 
electrostatic  energies  lower  than  the  specified  safe  level. 
Therefore,  the  net  charge  on  the  aircraft  is  kept  within 
acceptable  values. 

One  objective  of  this  report  is  to  provide  the 
appropriate  theoretical  analysis  of  dischargers  of  this 
type  and  of  its  several  components  in  order  to  facilitate 
and  improve  the  design  of  future  units.  Specifically, 
this  theoretical  information  is  intended  to  make  possible 
the  following  improvements  in  the  discharger: 

l.  Greater  corona  discharge  current 
at  a  given  voltage,  that  is,  lower 
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operating  voltage  and,  consequently, 
reduced  weight  and  complexity. 

2.  Lower  maintenance  requii  •‘ments. 

3.  Longer  component  life. 

Extensive  effort  has  been  directed  toward  finding 
a  method  for  the  reduction  of  operating  voltage  in  the 
discharger  system.  The  increase  of  current  dissipating 
capability  of  the  corona  electrode  would  permit  this 
reduction.  Consequently,  the  physics  of  corona  and 
electrostatic  fields  are  an  important  part  of  the  report. 
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SECTION  B:  CORONA  POINT  -  AIRCRAFT  FLSELAGE  CONFIGURATION 


Very  high  voltage  generators  are  required  to 
generate  the  corona  discharging  currents  essential  to 
maintain  a  low  helicopter  energy  level  under  charging 
conditions.  Recent  progress  in  the  state  of  the  art  in 
semiconductor  high  voltage  rectifiers  and  other  components 
permits  the  design  and  construction  of  highly  efficient, 
low  weight,  and  highly  reliable  high  voltage  D.C.  genera¬ 
tors  for  this  purpose. 

However,  in  addition  to  the  generation  of  a  high 
potential  on  a  corona  probe,  several  other  parameters 
must  be  considered  in  order  to  neutralize  an  aircraft 
electrostatic  charge  buildup.  Two  such  parameters  are: 

1.  Air  speed  vector  field  in  the  vicinity  of 
the  corona  probe. 

2.  Geometry  of  the  electrostatic  field 
generated  by  the  probe  potential. 

The  combination  of  these  two  factors  together 
with  the  probe  potential,  determine  the  discharging 
capabilities  of  the  probe. 

The  airspeed  vector  field  in  the  vicinity  of  the 
corona  probe  depends  on  the  location  of  the  probe  relative 
to  the  aircraft  fuselage.  To  effect  minimum  recirculation, 
the  airspeed  vectors  should  point  away  from  the  fuselage 
as  illustrated  in  Figure  1  for  both  single  and  tandem 
rotor  helicopters. 

For  practical  reasons,  the  corona  point  tip  should 
point  away  from  the  aircraft,  although  the  direction  of  the 
tip  is  not  expected  to  significantly  affect  the  direction 
of  motion  of  the  ions  emitted  from  the  probe  point. 

The  effectiveness  of  multiple  probes  connected  to 
the  same  high  voltage  generators  was  investigated  during 
flight  tests  performed  under  Contracts  DA  44-177-AMC-3(T) 
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Figure  1.  Typical  Corona  Probe  Configurations  for 
Single  and  Tandem  Rotor  Helicopters. 


and  DA  44-177-AMC-114(T) .  During  the  tests  conducted 
under  Contract  DA  44-177-AMC-3(T) ,  two  probes  which  were 
connected  to  the  same  high  voltage  generators  were  placed 
on  the  opposite  sides  of  the  aircraft  and  as  such  were 
electrically  shielded  from  each  other  by  the  aircraft 
structure.  It  was  found  that  for  the  same  generator 
voltage,  the  corona  current  of  both  probes  was  about 
double  that  of  a  single  probe.  On  the  other  hand,  under 
Contract  DA  44-177-AMC-114(T) ,  the  two  probes  were  mounted 
close  together  at  the  extremities  of  an  externally  mounted 
high  voltage  generator.  The  corona  current  for  this  case 
was  found  to  be  essentially  the  same  as  chac  of  a  single 
probe  at  the  same  voltage.  It  may  be  concluded,  therefore, 
that  substantial  separation  and  shielding  is  required  to 
obtain  increased  current  by  use  of  multiple  probes.  This, 
however,  implies  increased  system  weight  and  larger  air 
drag,  and  considering  the  present  state  of  the  art,  it 
is  much  more  efficient  to  increase  the  voltage  required  to 
obtain  the  required  discharge  current. 

The  determination  of  the  electrostatic  field 
distribution  near  an  active  corona  probe  installed  on  an 
aircraft  presents  a  classical  problem  of  potential  field 
theory.  A  generalized  solution  to  the  problem  will  be  of 
little  value,  because  the  inclusion  of  boundary  conditions 
closely  representing  aircraft  fuselage  configurations  will 
make  the  solution  to  the  problem  extremely  difficult.  It 
is  possible,  however,  to  find  a  solution  to  the  problem 
in  some  simplified  geometries,  such  as  the  installation 
of  a  corona  probe  on  an  infinitely  flat  surface.  This 
configuration  obviously  is  not  the  one  actually  encountered 
in  aircraft.  However,  the  field  distribution  near  the 
corona  probe  tip,  will  not  be  very  different  from  the 
field  distribution  existent  in  an  actual  installation.  The 
basis  for  this  statement  is  the  fact  that  near  the  corona 
tip,  the  main  field  component  is  due  to  the  tip  charge, 
the  surface  charge  on  the  airframe  introducing  only  a 
minor  effect  on  the  field  in  the  tip  area. 
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The  computation  of  the  electrostatic  field  of  a 
probe  placed  in  front  of  an  infinitely  conducting  plane 
can  be  achieved  by  assuming  a  charge  of  equal  magnitude 
and  opposite  polarity  located  at  a  point,  P2>  symmetrical 
to  the  corona  probe,  P^,  with  respect  to  the  infinite 
plane.  The  locations  of  the  points  P^  and  Po  are  shown  in 
Figure  2.  The  twin  charge  configuration  renders  an  electro¬ 
static  field  distribution  in  the  probe  side  of  the  plane, 
identical  to  the  field  distribution  that  would  result  by 
the  application  between  the  probe  and  the  plane  of  a 
potential  of  one-naif  that  assumed  between  the  points 
P^  and  P9. 


To  prove  this  similarity,  let  us  assume  there  is 
no  conductive  plane  between  the  charges  at  P^  and  P?.  The 
field  vector  at  point  "O”  is  E  and  is  parallel  to  the  line 
between  the  charge  points.  At  any  other  point  B  on  the 
axis  of  symmetry  the  field  vector  Eg  is  parallel  to  the 
line  between  P^  and  P2  because  its  two  components  E»i  and 
Eg2  have  the  same  magnitude  and  are  symmetrical  witn 
respect  to  a  line  parallel  to  P^  P2  through  point  B.  Hence, 
the  plane  of  symmetry  is  an  equipotential  surface  in  the 
field  created  by  the  charges  at  P^  and  P2,  because  the 
motion  of  a  charge  on  it  can  be  achieved  without  giving 
or  receiving  work  on  account  of  the  perpendicularity  be¬ 
tween  any  charge  path  within  the  plane  and  the  field  vector 
direction  on  it.  Thus,  the  resultant  field  on  the  probe 
side  of  this  conductive  plane  will  be  identical  to  the 
field  distribution  created  by  the  two  charges  at  points 
P^  and  P2. 


The  field  at  a  point  A,  whose  coordinates  are 
R  and  h  can  be  now  computed  as  follows. 


The  field  due  to  the  charge  at  P-.  is 

a  % 


4  s 


y/re. 


Pj_ 


¥n  Vi 


1 


but 


z  A  ft  1  •  (fi  "A)  JTj 


(B-l) 


(B-2) 


7 


Figure  2. 
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where  and  riy  are  two  unit  vectors,  n,  being  perpendi¬ 
cular  to  line  P1P2  and  being  parallel  to  line  Pj^. 


The  field  due  to  the  charge  at  ?2  *s 
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and  the  total  field  at-  a  will  be 
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(B-7) 


The  two  coefficients  of  equation  (B-7)  are  the 
magnitude  of  the  two  components  of  the  electrostatic  field 
parallel  and  normal  to  the  conductive  plane.  The  symmetry 
around  the  corona  probe  axis  permits  the  computation  of 
these  two  field  vector  components  at  any  point  in  the 
space  between  the  probe  and  the  plane. 
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The  field  E  in  equation  (B-7)  is  presented  as 
a  function  of  the  probe  charge  qD  . 


The  relationship  between  the  probe  charge  q0  and 
its  potential  with  respect  to  the  airframe  can  be  established 
by  using  again  the  symmetrical  charge  analogy.  By  symmetry, 
it  is  obvious  that  the  potential  between  both  charges  is 
double  that  of  the  potential  between  one  charge  and  the  plane 
of  symmetry.  The  potential  between  both  charges  can  be  com¬ 
puted,  as  the  work  required  tc  move  a  unit  of  charge  between 
them  against  the  electrostatic  field  of  force.  This  potential 
is  easily  obtained  for  R  =  0, 


So 


y/r  <£>•  A* 


HA 
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where  Vp  is  the  potential  between  the  charges, 


(B-8) 


qQ  is  the  charge  on  the  probe  tip,  and 


R^.  is  the  radius  of  the  probe  tip. 


The  force  acting  on  a  particle  with  charge  q  is 
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SECTION  C:  CORONA  DISCHARGE  IN  STILL  AIR 


An  account  of  the  processes  occurring  in  corona 
discharge  is  presented  as  a  basis  for  the  mathematical 
analysis  of  corona  given  in  the  next  section.  This 
analysis  provides  a  basis  for  the  design  of  efficient 
electrostatic  discharging  systems.  Accordingly,  onset 
conditions,  characteristics,  and  comparisons  of  positive 
and  negative  corona  are  given  in  this  section.  This  infor¬ 
mation  was  derived  from  a  study  of  ion  generation  in 
atmospheric  corona  discharge  (References  6,  11,  12  and  18). 

Corona  is  the  name  given  to  electrical  conduction 
in  a  gas  just  before  arc-over.  It  is  distinct  from  other 
stages  of  conduction  in  that  avalanche  ionization  takes 
place,  and  the  discharge  is  self-sustaining,  that  is,  no 
external  ionization  source  such  as  cosmic  radiation  is 
necessary  for  its  continuance.  Whatever  the  type  of  gas 
discharge,  ionization  is  always  present  since  without  it 
a  gas  is  nonconductive.  A  small  amount  of  natural 
ionization  due  to  cosmic  rays  and  radioactive  trace 
materials  suffices  to  trigger  a  discharge  when  a  voltage 
is  applied.  An  important  feature  of  gaseous  conduction 
is  the  nature  of  the  deionization  process,  that  is,  the 
loss  of  ions  constituting  the  discharge.  This  loss  can 
occur  in  three  ways:  recombination  of  electrons  and 
positive  ions  to  neutral  atoms,  attachment  of  electrons 
to  immobile  molecules,  and  diffusion  out  of  the  discharge 
region. 

ONSET  CONDITIONS  AND  CHARACTERISTICS  OF  NEGATIVE  CORONA 


There  are  three  stages  of  discharge  in  the 
atmosphere  which  occur  in  sequence  as  the  negative  high 
voltage  applied  to  the  discharge  electrode  is  increased 
to  the  corona  threshold  value. 

Townsend  Discharge 

The  values  of  negative  voltage,  being  first  very 
small,  allow  recombination  to  occur  in  the  discharge. 

This  causes  the  relationship  between  the  discharger  system 
electrode  voltage  and  discharge  current  to  be  initially 
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linear.  When  the  negative  voltage  is  further  increased,  the 
ions  of  the  discharge  are  drawn  to  the  electrode  before  any 
recombination  can  occur.  This  causes  current  to  remain  con¬ 
stant  over  a  small  range  of  voltage.  As  voltage  is  still 
further  increased,  the  corona  discharge  current  begins  to  in¬ 
crease  sharply.  This  large  increase  occurs  because  the  ex¬ 
ternally  caused  ions,  that  is,  ions  resulting  from  cosmic 
rays  or  radioactive  trace  materials,  have  acquired  sufficient 
energy  to  cause  further  ionization.  The  discharge,  however, 
is  not  yet  self-sustaining. 

Electrical  Breakdown 


The  initial  phase  of  breakdown  is  ionization  by 
the  process  of  externally  caused  ionization.  This  con¬ 
dition  is  termed  ,?electron  avalanche".  The  positive  ions 
resulting  from  the  avalanche  generate  new  electrons  by 
ionizing  air  molecules  while  being  drawn  to  the  negative 
electrode.  Also,  electrons  are  generated  by  secondary 
emission  when  these  positive  ions  strike  the  negative 
electrode.  A  condition  is  reached  as  voltage  is  increased 
in  which,  on  the  average,  the  products  of  one  avalanche 
produce  at  least  one  electron  capable  of  starting  a  second 
avalanche  as  large  as  the  first.  An  equilibrium  between 
ionization  and  deionization  follows  and  the  discharge,  being 
independent  of  external  sources  of  ionization,  becomes  self- 
maintaining.  This  condition  is  termed  "electrical  break¬ 
down"  and  is  present  in  corona  discharge 

The  ability  of  an  electron  to  cause  ionization  and 
subsequent  breakdown,  and  therefore  contribute  to 
the  discharge  current,  depends  on  the  energy  it  receives 
from  the  electric  field  between  collisions  with  the  air 
molecules.  This  energy  is  given  by 

£  *  £\e  (c-i) 

=  electric  field  strength, 

=  average  distance  the  electron  travels 
between  collisions, 

=  ion  kinetic  energy  increase,  and 
=  electrical  charge  of  an  electron. 


where 


E 

A 
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The  quantity  A  is  directly  proportional  to  tem¬ 
perature  and  is  inversely  proportional  to  pressure.  Hence, 
at  constant  temperature 

£  *  A,  £2-  (C-2) 

where  A,  =  mean  free  path  of  an  electron  at 

1  mm  Hg  and 

p  =  atmospheric  pressure  in  units  of  mm  Hg. 

Therefore,  the  electrical  quantity  of  importance  in  corona 
discharge  current  calculations  at  constant  temperature  is 
E/p. 


The  expression  for  discharge  current  density 
provides  the  mathematical  condition  for  corona  onset. 

This  expression  can  be  calculated  by  the  use  of  Townsend *s 
first  and  second  ionization  coefficients  (Reference  19). 

They  are,  respectively, 

a  =  number  of  ionizing  collisions  per  unit 
length  of  path  in  the  direction  of  the 
field,  and 

Y  -  number  of  new  electrons  produced  per 
positive  ion  due  to  any  or  all  of  the 
mechanisms  of  ionization  by  positive 
ions,  secondary  emission  at  the  nega¬ 
tive  point,  and  photo -ionization. 

The  quantity  Of  depends  on  the  energy  of  the  colliding 
electron,  that  is  E/p,  and  the  number  of  collisions  per 
unit  path  length  which  is  directly  proportional  to  pressure, 
that  is, 

cx.  -  (c-3) 

Therefore,  °<7p  depends  on  E/p  only.  The  second  ioniza¬ 
tion  coefficient  f  depends  on  the  material  of  the  point 
and  the  gas.  It  takes  into  account  the  emission  of  se¬ 
condary  electrons  by  the  positive  ions  which  is  not  signi¬ 
ficant  until  the  latter  stages  of  breakdown. 
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Before  calculating  current  density,  the  electrode 
geometry  must  be  specified.  The  discharger  system  electrode 
geometry  is  most  closely  simulated  by  a  point  conductor 
(corresponding  to  the  corona  electrode)  suspended  over  an 
infinite  plane  (corresponding  to  the  aircraft)  at  a  speci¬ 
fied  distance.  To  a  close  approximation,  the  field  varies 
only  with  the  position  between  point  and  plane.  The 
distance  from  the  point  conductor  along  the  perpendicular 
to  the  infinite  plane  from  the  point  is  the  X  coordinate. 

Let  the  rate  of  natural  or  artificial  ionization 
triggering  the  breakdown  be  distributed  uniformly  in  the 
inter-electrode  space  with  a  magnitude  of  ions  per 
unit  volume  per  second.  In  traversing  an  elementary 
volume  of  unit  area  and  length  dx,  the  colliding  electrons 
produce  an  amount  dn  of  new  ions  per  second  which  is 
proportional  to  n^,  the  number  of  electrons  entering  this 
volume  per  second,  and  the  length  dx.  The  proportionality 
constant  is  Townsend's  first  ionization  coefficient  o< 

Thus,  if  recombination  and  diffusion  are  neglected  and  if 
°<  is  assumed  constant  throughout  the  region  under  con¬ 
sideration 


</a  *  (*n,  *  9<)  o(> c 


(C-4) 


If  nQ^  is  the  number  of  electrons  emitted  at  the  negative 
point  (at  x  =  0)  by  cosmic  or  artificial  irradiation, 
integration  of  equation  (C-4)  gives 


H  X 

n,  (x)  *  n0/  e  ♦ 


-*-(e 

oC  \ 


*  ') 


(C-5) 


Since  oC &  10  for  E/p  =  60  volts/centimeter  at  760  mm 
Hg  (Reference  6,  page  150),  there  follows  that  for  X  »0 

1.  Hence  equation  (C-5)  can  be  reduced  to 


the  term  e  » 


n,(x) 


m  — 

n  ♦  ^ 
•c  _ 


(C-6) 


Equation  (C-5)  is  the  electron  current  density  in  the 
units  of  electrons  per  second  per  square  centimeter  at  an 
arbitrary  position  in  the  gaseous  discharge.  No  considera¬ 
tion  of  secondary  emission  of  electrons  is  made  in  ob¬ 
taining  this  expression.  In  order  to  obtain  the  con¬ 
dition  for  corona  onset,  this  expression  must  be  modified 
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to  give  the  electron  current  density  when  secondary 
emission  is  considered.  This  modification  involves  re¬ 
placing  the  quantity  nQ^  in  equation  (05)  by  its 
equivalent,  nc,  when  secondary  emission  is  present 
(Reference  6).  Let 


na  =  number  of  electrons  reaching  anode  per 
unit  area  per  second, 


no2 


n 


c 


c 


number  of  electrons  per  unit  area  per 
second  emitted  from  negative  point  due 
to  cosmic  rays,  et  cetera, 

ions  per  unit  volume  per  second  formed 
in  inter-electrode  space  due  to  cosmic 
rays,  et  cetera, 

total  number  of  electrons  per  unit  area 
per  second  emitted  from  negative  electrode 
due  to  cosmic  rays,  et  cetera  and  to  secon¬ 
dary  emission,  and 

distance  from  negative  electrode  to 
anode . 


Then,  na  -  nc  =  number  of  positive  ions  formed  in  the 

atmosphere  and 


r\c  '  rioz  *  f  (”*  "  He) 


(C-7) 


nc  *  -A.*..*  rn* 

/ 1  r 


(C-8) 


The  quantity  nc  is  the  equivalent  of  n0^  when  secondary 
emission  is  considered  and  replaces  it  in  the  expression 
for  electron  current  density  at  any  position  x  in  the 
discharge  equation  (C-6).  Therefore,  with  secondary 
emission  considered,  the  electron  current  density  at  the 
anode  is 

/?.  *  («c  *  e"c  (c-9) 
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(C-10) 


Replacing  n_  by  using  equation  (C-8)  gives 

-  \zaJHHa.  ♦  Ji]  e*e 
/  <* 

Removing  n  from  the  right  hand  member  of  equation  (C-10) 
we  have  A 

n«  •  K  *•£■('*')}  (c'u) 

/  -  +  r 

The  total  current  density  in  gaseous  discharge  equals  the  sum 
of  the  positive  ion  and  electron  current  densities.  The  law 
of  current  continuity  states  that  total  current  density  in 
any  gaseous  discharge  is  constant  for  any  position  in  the 
discharge  when  steady  state  conditions  exist.  Therefore, 
the  total  current  density  equals  the  electron  current  density 
na  (multiplied  by  e,  the  electron  charge)  at  the  position 
in  the  discharge  where  the  positive  ion  current  is  zero. 

This  point  is  the  anode.  Therefore, 

Jt  .  nKe  •  g  *  ~%c(‘  *_  *j\  e“c  (c-12) 

r  e*c 

where  Jt  =  total  current  density  before  breakdown  and 
e  =  electronic  charge. 

The  quantity  has  been  neglected  in  comparison  to 

in  the  denominator  of  equation  (C-ll) .  Equation  (C-12)  is 
of  practical  use  in  predicting  magnitudes  of  pre-breakdown 
discharge  current. 
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It  indicates  that  selecting  the  electrode  material  for 
proper  y  can  increase  discharge  current  from  the  electrode. 

The  denominator  of  equation  (C-12)  gives  the  con¬ 
dition  for  negative  corona  onset.  The  condition  is  that 
the  current  density  Jt  increase  without  limit.  This  occurs 
if  the  denominator  of  equation  (C-12)  goes  to  zero,  that 
is  if 

w  -  C 

6  e  *  1  (C-13) 


Equation  (C-13)  implies  that  corona  commences  at 
that  value  of  electrical  field  which  causes  at  least  one 
new  electron  to  be  emitted  from  the  negative  electrode 
when  it  is  bombarded  by  the  positive  ions  present  in  an 
avalanche.  The  discharge  is  termed  ’’self-sustaining”  at 
this  stage.  The  exact  value  of  electric  field  at  which 
corona  occurs  in  air  at  constant  temperature  and  pressure 
is  a  function  of  corona  electrode  radius  of  curvature.  A 
low  corona  threshold  value  is  very  favorable  in  that  dis¬ 
charge  current  is  then  larger  at  a  given  electrode  voltage 
as  shown  by  equation  (D-50).  Reference  19,  page  153,  gives 
data  on  the  corona  onset  values  of  electric  field  for  nega¬ 
tive  points  with  radii  ranging  from  0.00023  inches  to  0.018 
inches  at  different  point  to  plane  spacings  (gap  lengths) . 

By  the  term  ’’plane”  we  make  reference  to  the  infinite  plane 
approximation  to  the  aircraft  fuselage  discussed  in  Section 
B.  Thus,  the  aircraft  fuselage  is  considered  to 
have  the  same  effect  ou  the  corona  current  of  the  electro¬ 
static  discharger  system  as  the  plane  used  in  the  labora¬ 
tory  experiments  described  in  Reference  19.  The  data  show 
that,  for  a  given  point  radius  of  curvature,  the  corona 
onset  electric  field  has  become  unchanging  with  respect 
to  change  in  gap  length  for  gap  lengths  greater  than  approxi¬ 
mately  3.14  inches  (8  centimeters).  Also,  for  a  given  gap 
length,  the  value  of  corona  onset  electric  field  decreases 
as  radius  of  curvature  increases  approaching  a  constant 
value,  Ep  ~  62  kilovclts/centimeter  at  a  point  radius  of 
0.038  inch  for  large  gap  lengths.  The  values  of  corona 
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onset  electric  field  in  this  data  were  found  by  using  the 
applied  voltage  and  appropriate  physical  dimensions  in 
the  laboratory  experiment  in  t.'ie  expression  for  electric 
field  at  the  corona  point  for  hyperboloidal  point  and 
opposing  plane  geometry  given  by 


where  c 

Vc 


F  .  j  _  (C-14) 

*  (-%) 

distance  between  corona  point  and  plane, 

corona  probe  threshold  voltage, 

radius  of  curvature  of  the  point,  and 

field  intensity  at  the  discharging 
electrode, 


The  preceding  expression  and  data  can  be  used 
to  give  data  meaningful  to  the  electrostatic  discharger 
system.  The  basis  of  applicability  is  that  as  discussed 
above,  for  gap  lengths  greater  than  3.14  inches  (8  centi¬ 
meters)  the  value  of  corona  onset  electric  field  is  con¬ 
stant  for  a  given  corona  point  radius.  Thus,  the  potential 
Vc  with  respect  to  the  aircraft  necessary  for  onset  of  the 
corona  current  in  the  electrostatic  discharger  system  can 
be  calculated  from  the  relationship  (C-14)  using  the  value 
of  point  to  plane  (fuselage)  spacings,  c  (greater  than  3.14 
inches),  the  corona  point  radius  of  curvature  R^,  and  the 
value  of  onset  electric  field  appropriate  to  gap  length 
(given  in  the  data  from  Reference  19),  Figure  3  shows 
curves  resulting  from  such  computations  for  Vc.  These 
curves  give  the  relationship  between  negative  corona 
threshold  voltage  and  electrode  tip  radius  of  curvature 
for  three  different  distances  between  aircraft  and  corona 
electrode.  The  symbols  used  in  these  curves  (squares,  et 
cetera)  indicate  the  values  of  point  radius  of  curvature  and 
point -fuselage  spacing  for  which  values  of  Vc  were  actually 
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Negative  Corona  Threshold  Voltage,  Vc,  Volts 


Radius  of  Curvature,  Rt,  Inches 


Figure  3.  Minimum  Negative  Probe  Voltage  Required 
for  Initiating  Corona  as  a  Function  of 
Electrode  Tip  Radius  and  Spacing  from 
Aircraft. 
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calculated.  The  four  maximum  and  minimum  points  on  the 
curves  result  from  the  fact  that  field  intensity  Ep 
varies  nonlinearly  both  with  Rt  and  c,  as  indicated  in 
Reference  19,  page  153.  The  curves  indicate  that  the 
corona  threshold  voltage  is  approximately  independent 
of  corona  electrode  tip  radius  for  values  of  radius 
between  0.002  and  0.006  inch.  Corona  point  deterio¬ 
ration  because  of  positive  ion  bombardment  (sputtering) 
makes  it  impossible  to  maintain  a  tip  radius  of  less 
than  0.001  inch  in  actual  operation  (Reference  18). 

Thus,  referring  to  Figure  3,  the  best  corona  electrbde 
tip  radius  of  curvature  is  between  0.002  and  0.006  inch, 
and  is  therefore  a  noncritical  dimension. 

Corona  Discharge 

The  term  "corona  discharge”  actually  signi¬ 
fies  glow  discharge  on  a  curved  electrode  (Reference  9). 

Thus,  corona  physics  is  that  of  the  familiar  glow  discharge. 
This  type  of  gaseous  conduction  is  characterized  by  low 
voltage  and  is  independent  of  current  density.  Another 
characteristic  of  the  corona  discharge  is  a  voltage  drop 
concentrated  over  a  short  distance  from  the  cathode  about 
equal  in  magnitude  to  the  breakdown  potential  of  air. 

'"his  concentration  of  voltage  drop  is  caused  by  the 
tendency  of  the  positive  ions  to  drift  toward  the  cathode 
and  form  a  virtual  anode  there.  Microscopic  investi¬ 
gation  of  the  voltage  drop  region  shows  it  to  be  composed 
of  a  series  of  relatively  light  and  dark  regions  called 
the  pston  dark  space.  The  positive  column, a  region  of 
small  voltage  gradient  due  to  small  net  space  charge 
and  masked  cathode,  comprises  the  remainder  of  the 
discharge.  In  normal  operation  of  the  discharger  system, 
the  corona  discharge  is  confined  to  a  region  only  a  few 
electrode  tip  radii,  (Reference  18).  The  discharge  region 
extends  to  the  aircraft  only  if  arc-over  occurs. 

The  phenomena  occurring  in  corona  discharge  as 
the  negative  electrode  voltage  is  increased  from  the  corona 
threshold  value  are  described  as  follows.  The  rate  of  change 
of  potential  with  distance  ceases  to  be  a  constant  as 
the  effects  of  space  charge  caused  by  ionization  begin 
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to  set  in.  As  voltage  is  increased  further,  the  cathode 
fall  of  potential  becomes  a  constant  determined  by  the 
electrode  material  (approximately  300  volts),  and  the 
current  emitting  area  of  the  electrode  increases  causing 
discharge  current  to  increase.  The  current  density, 
however,  remains  constant  at  this  stage.  The  remainder 
of  the  corona  electrode  high  voltage  developed  by  the 
system  is  dropped  from  the  periphery  of  the  corona  dis¬ 
charge  to  the  aircraft  surface.  With  further  increase 
of  negative  voltage,  the  current  density  begins  to  increase 
and  arcs  or  streamers  occur  between  the  electrode  and  the 
virtual  anode.  Destruction  of  the  discharge  electrode 
follows  because  of  positive  ion  bombardment  and  heavy 
current  loading  (Reference  18). 

The  ionization  processes  in  the  atmosphere  cause 
the  presence  of  electrons,  positive  oxygen  ions,  negative 
oxygen  ions,  positive  nitrogen  ions,  and  ionizing  photons 
as  constituents  of  the  discharge.  The  primary  processes 
are 


Oa  — 

♦  °1 

+  electron 

(C-15) 

Nx  — 

-  N* 

+  electron 

(C-16) 

In  closed  spaces,  the  nitrogen  gaseous  compounds 
are  noxious.  The  corona  current  in  the  atmosphere  has  a 
unique  pulsed  characteristic  (Reference  12)  because  of 
the  presence  of  an  electronegative  gas  (oxygen).  Neutral 
molecules  of  electronegative  gases  have  a  great  affinity 
to  attach  electrons  in  regions  of  small  E/p,  that  is, 
where  electrons  have  small  velocities.  The  resulting  heavy 
immobile  negative  ions  accumulate  and  eventually  choke-off 
the  corona  current  momentarily,  thereby  causing  the  pulsed 
characteristic.  Corona  current  in  pure  nitrogen  and  hydro¬ 
gen  gases  which  are  not  electronegative,  are  much  greater 
in  magnitude  and  are  not  pulsed.  These  gases,  once  the 
steady  state  is  reached,  have  space  charge  distributions 
depending  on  position  only.  For  the  Trichel  corona  current 
pulse,  named  after  its  discoverer,  the  space  charge  distri- 
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button  depends  on  both  position  relative  to  the  corona 
electrode  and  time.  The  characteristics  of  the  Trichel 
pulse  are  important  because  the  system  discharge  current 
is  the  time  average  of  the  pulsed  current. 

The  steps  in  the  negative  pulse  formation  are  as 
follows.  A  random  positive  ion  due  to  cosmic  radiation 
is  accelerated  toward  the  negative  electrode.  Initially, 
the  field  is  high  enough  to  accelerate  the  ion  to 
sufficient  energy  to  cause  the  emission  of  a  secondary 
electron  when  it  reaches  the  electrode.  This  electron 
produces  a  dense  cloud  of  positive  ions  by  the  Townsend 
avalanche  mechanism  as  it  is  repelled  from  the  negative 
electrode.  This  slow  moving  positive  ion  space  charge 
reduces  the  field  of  the  negative  electrode  at  its  outer 
periphery  where  it  is  already  weak  due  to  distance  from 
the  electrode.  Consequently,  the  repelled  electrons  lose 
energy  rapidly  and  form  negative  ions  by  attachment  to 
neutral  oxygen  molecules.  The  accumulation  of  these  nega¬ 
tive  immobile  ions  reduces  the  field  which  had  been 
drawing  the  positive  ions  to  the  electrode.  It  is  this 
effect  which  tends  to  choke  off  the  current. 

The  decelerated  positive  ions  proceed  without 
producing  further  ionization  to  the  electrode  where  they 
form  the  current  pulse  and  contribute  to  the  discharge 
of  the  negatively  charged  aircraft.  By  the  time  the 
last  few  positive  ions  approach  the  electrode,  the 
negative  heavy  ion  space  charge  may  have  been  cleared 
sufficiently  by  the  viscous  downwash  forces,  that  the 
field  is  high  enough  to  cause  another  repetition  of  the 
pulse.  It  Is  unlikely  that  free  electrons  attach  to 
form  ions  in  regions  of  the  discharge  where  they  have 
an  average  energy  corresponding  to  that  causing  ex¬ 
citation  and  ionization,  that  is,  where  E/p  exceeds  90 
volts  per  centimeter  per  millimeter  Hg,  since  detach¬ 
ment  of  electrons  from  the  neutral  oxygen  molecules 
occurs  at  these  values  of  E/p. 
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The  corona  electrode  performance  is  improved  due 
to  its  location  in  the  rotor  downwash  region.  The  airstream 
tends  to  prevent  accumulation  of  the  heavy  negative  ions 
around  the  electrode.  The  corona  current  is  therefore 
greater  since  less  choke-off  occurs. 

Parameters  of  the  negative  corona  current  pulse 
such  as  pulse  repetition  rate,  duration,  and  rise  time 
depend  on  the  following: 

1.  Electrode  configuration 

2.  Electrode-aircraft  separation 

3.  Material  of  the  electrode 

4.  Environmental  conditions 

5.  Applied  voltage 

The  bases  for  the  above  statement  are  the  following  consid¬ 
erations  which  also  indicate  how  these  conditions  affect 
the  corona  current  pulse  parameters.  By  definition  the 
average  corona  current  is  the  net  amount  of  charge  en¬ 
tering  the  corona  point  divided  by  the  time  interval 
required  for  this  charge  to  enter  the  point.  In  the 
pulsed  nature  of  negative  corona  current  as  qualitatively 
described  in  the  preceding  paragraph,  a  fixed  amount  of 
charge  enters  the  point  in  a  specific  small  time  interval 
(pulse  duration  time).  (In  general,  for  voltages  near 
corona  onset,  the  Trichel  pulses  are  somewhat  random  and 
irregular.  As  the  electrode  voltage  is  increased,  the 
pulses  become  regular.)  Thus,  in  mathematical  terras,  the 
average  corona  current,  I,  Is  given  (Reference  16)  by 


i  •  Ft  em 


(C-17) 


where 


I  is  corona  discharge  current, 

Fr  is  pulse  repetition  rate  (reciprocal  of 
pulse  duration  time) , 
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m  is  number  of  ions  involved  in  one  current 
pulse,  and 

e  is  electronic  charge. 

The  quantity  m  is  given  in  the  Townsend  avalanche 
mechanism  theory  (Reference  19)  as 


m 


(C-18) 


where  nQ  is  the  number  of  electrons  emitted  at 

the  negative  corona  electrode, 

a  is  Townsend’s  first  ionization  coefficient,  and 

c  is  distance  from  negative  corona  electrode 
to  virtual  anode. 

The  integral,  dx,  is  not  exactly  evaluable 
because  of  the  complicated  dependence  of  on  x  through 
the  relationship  of  E  with  x  (References  19  and  11),  and 
recalling  that  *>c  -  p/(E/p),  see  equation  (C-3).  The 
latter  relationship  differs  from  the  basic  E  versus  distance 
relationship  given  in  Section  B  because  of  the  effects  of 
space  charge  due  to  atmospheric  ionization.  The  exact  re¬ 
lationship  can  be  obtained  through  the  differential  equa¬ 
tions  of  Section  D,  page  56.  The  literature  (Reference 
11)  describes  the  exact  equations  as  solvable  by  computer 
(differential  analyzer).  In  practice,  the  number  of  ions 
involved  in  a  pulse  is  obtained  through  graphical  inte¬ 
gration  of  the  current  pulse  shape  to  determine  the  in¬ 
cluded  area  which  equals  the  amount  of  charge  (number  of 
ions  multiplied  by  electronic  charge,  e). 

The  faster  the  ions  involved  in  the  pulse  move, 
the  shorter  will  be  the  time  duration  of  the  pulse  and 
the  greater  the  repetition  rate,  Fr,  will  become.  Ion 
velocity  is  determined  by  electrical  field  intensity,  E, 
being  directly  proportional  to  it  with  constant  of  pro¬ 
portionality,  kj,  the  ion  mobility  (discussed  in  Section 
D,  page  51).  Thus,  the  following  equation  can  be  written 
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s  kr  v  •  £ 


(C-19) 


where 


v  is  ion  velocity, 

Fr  is  pulse  repetition  rate, 

kr  is  the  constant  of  proportionality  between 

Fr  and  v, 

is  the  ion  mobility  constant  of  proportionali¬ 
ty  ,  and 

E  is  electric  field  intensity. 

Hence,  equation  (C-17)  can  be  rewritten  as  follows: 


f  C/>f(£4>')ctx  (C-20) 

x  «  AnM  eJ»  ' 

where  A.  is  current  emissive  surface  ..rea  of  corona 

electrode , 

noA  is  number  of  ions  emitted  at  corona  electrode 
per  unit  area,  and 

p,  ky,  k^,  E,  e,  and  x  have  been  defined  previously. 


Equation  (C-20)  shows  quantitatively  the  way  in 
which  the  five  conditions  listed  on  page  23  influence 
corona  current.  Electrode  configuration  determines  A  and 
the  value  of  E  at  any  point  in  the  discharge.  Electrode 
aircraft  separation  determines  the  value  of  E  at  points 
in  the  discharge  in  accordance  with  the  considerations  in 
Section  B.  Material  of  the  electrode  affects  n^  through 
the  mechanism  of  secondary  emission  from  the  electrode. 
Also,  the  material  of  the  electrode  affects  its  configura¬ 
tion  in  consequence  of  durability  considerations.  En¬ 
vironmental  conditions  determine  p  and,  through  tempera¬ 
ture,  affect  ion  mobility.  Applied  voltage  partially  de¬ 
termines  E  (refer  to  equations  on  page  56) .  Since  an  in¬ 
crease  of  the  applied  potential  accelerates  the  Trichel 
pulse  cycle,  the  pulse  repetition  rate,  Fr,  is  increased 
with  an  increase  of  applied  voltage.  The  number  of  ions 
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involved  in  one  current  pulse,  m,  remains  constant  as  the 
negative  electrode  voltage  is  increased  from  the  corona 
threshold  value.  Hence,  near  onset,  negative  corona  current 
is  directly  proportional  to  applied  voltage.  This  holds 
over  a  considerable  range  in  voltage  but  eventually,  as 
electrode  voltage  is  increased,  there  is  a  departure  from 
the  linear  current.  This  indicates  a  decrease  in  the  number 
of  ions,  m,  per  pulse.  A  Trichel  pulse  repetition  rate  of 
10  megacycles /second  has  been  observed  before  the  applied 
voltage  reached  the  arc-over  point  and  100  megacycles/second 
is  theoretically  possible  (Reference  19) .  The  number  of  ions 
m  for  a  point  of  0.2-mm  radius  is  about  10^  ions.  Thus, 
application  of  equation  (C-17)  indicates  that  a  current  of 
1  microamp  corresponds  to  a  pulse  repetition  rate  of  6 
kilocycles  per  second.  A  current  of  about  50  microamps 
can  be  obtained  for  an  applied  voltage  of  4U  kilovolts  in 
a  point-plane  arrangement  with  a  gap  of  about  6  inches. 

The  near  proximity  of  the  plane  to  the  corona  electrode 
increases  the  corona  current  significantly  because  of  the 
resulting  increased  field  strength  at  the  electrode  as 
explained  in  Section  D.  Of  course,  current  going  to  the 
infinite  plane  is  analogous  to  discharge  current  returning 
to  the  aircraft.  An  actual  discharger  system  permitting 
this  to  occur  would  not  be  acceptable.  The  action  of  the 
air stream  tends  to  prevent  this  by  diverting  this  current 
from  the  aircraft,  if  the  airspeed  is  sufficiently  high, 
as,  for  example,  for  a  very  high  speed  aircraft.  For 
helicopters,  however,  the  maximum  available  airspeed  is  not 
sufficient  to  permit  the  use  of  this  approach  for  the 
design  of  the  corona  probe.  Reference  14  gives  values  of 
rise  time  and  duration  for  the  negative  Trichel  pulse. 

Typical  values  at  atmospheric  pressure  are  risetimes  of 
12  nanoseconds  and  pulse  duration  of  30  nanoseconds  for  a 
0.015-inch  radius  electrode  tip  at  30  kilovolts. 


POSITIVE  CORONA  CHARACTERISTICS 


Two  types  of  positive  corona,  depending  on  at¬ 
mospheric  characteristics  and  applied  voltage,  may  occur. 
These  are  burst  pulse  corona  similar  to  the  Trichel  pulse, 
and  streamer  (or  brush)  corona  (Reference  19).  They  differ 
in  appearance  markedly.  By  considerations  similar  to  those 
used  in  finding  the  negative  corona  onset  criterion,  the 
condition  for  intermittent  burst  pulse  corona  may  be  de¬ 
rived  as 
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(C-21) 


r  / 


The  significance  of  quantities  in  this  expression  are  dis¬ 
cussed  below.  Here  as  in  the  negative  corona  case, 


e 


is  the  number  of  ions  in  an  electron  avalanche.  In  the 
case  of  positive  corona,  the  avalanche  is  due  to  ioniza¬ 
tion  by  a  naturally  occurring  trigger  electron  as  it  is 
accelerated  toward  the  corona  electrode.  The  ions  which 
have  provided  the  electrons  for  the  avalanche  constitute 
a  positive  space  charge  when  the  avalanche  has  entered  the 
positive  electrode.  Similar  to  the  Trichel  pulse  pheno¬ 
menon  of  negative  corona  discharge,  this  positive  space 
charge  causes  burst  pulses  of  corona  discharge  current. 

The  burst  pulse  duration  can  be  from  tens  to  hundreds  of 
microseconds.  In  contrast  to  the  negative  corona  mechanism, 
the  appearance  of  new  electrons  is  not  due  to  secondary 
emission,  but  due  to  photoionization.  This  occurs  during 
the  last  portion  of  the  path  traveled  by  the  avalanche 
electrons  which,  before  entering  the  positive  electrode, 
create  a  number  of  photons  of  high  energy  ultraviolet 
light  by  impacts  with  atoms.  The  number  of  photons  pro¬ 
duced  is  proportional  to 

e  J° 

and  is  therefore  given  by 


c 

otx 
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in  which  fp  is  a  constant  of  proportionality.  A  fraction, 

Q  ,  of  these  photons  escape  absorption  and  create  elec¬ 
trons  in  the  gas  within  a  distance  c  from  the  electrode 
by  photoionization.  Any  of  these  photoelectrons  can  cause 
a  new  avalanche  which  will  perpetuate  the  process.  The 
onset  criterion  in  equation  (021)  signifies  the  pro¬ 
duction  of  at  least  one  new  electron  to  start  a  new 
avalanche  and  therefore  self -maintenance  of  the  discharge. 

Streamers  occur  when  the  positive  ion  charge 
c.ensity  is  large  enough  to  act  as  an  extension  of  the 
positive  electrode.  Under  these  conditions,  electrons 
are  drawn  to  the  positive  ion  charge  density.  They 
ionize  air  molecules  causing  more  positive  space  charge. 

Thus  the  streamer  extends  itself  considerably  into  the 
gap  as  a  luminous  conducting  filament.  Streamers  appear 
only  at  low  currents  or  at  high  currents  just  before 
breakdown.  Consequently,  the  current  over  most  of  the 
gap  between  the  point  and  pLane  consists  of  randomly 
spaced  burst  pulses. 

Figure  4  shows  curves  of  positive  corona  thresh¬ 
old  voltage  as  a  function  of  radius  of  curvature.  These 
curves  were  calculated  similarly  to  those  for  negative 
corona  threshold,  that  is,  the  point -plane  approximation  was 
used.  The  remarks  made  for  negative  corona  with  respect 
to  dimension  of  electrode  tip  also  apply  for  positive  corona. 

COMPARISON  BETWEEN  POSITIVE  AND  NEGATIVE  CORONA 


Consideration  of  the  positive  and  negative  corona 
onset  conditions  with  the  values  of  if  and  ft  fp  appropiiate 
to  air  (  Y  &  lO’^j^Pfp^  10"^)  applied  seems  to  indicate 
that  the  negative  Trichel  pulse  onset  voltage  is  lower 
than  the  positive  burst  pulse  onset  voltag?  (Reference  12). 
It  would  also  be  expected  that  the  negative  onset  voltage 
which  depends  on  X  would  be  sensitive  to  the  electrode 
material  and  condition  while  the  positive  onset  voltage 
would  depend  on  the  gas  only.  In  air,  however,  the  two 
thresholds  are  markedly  close  together.  The  thresholds 
are  about  the  same,  even  when  corona  electrodes  of  different 
metals  such  as  platinum,  copper,  and  zinc  are  used.  The 
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Positive  Corona  Threshold  Voltage,  Vc,  Volts 


Radius  of  Curvature,  Rt,  Inches 


Figure  4.  Minimum  Positive  Probe  Voltage  Required 
for  Initiating  Corona  as  a  Function  of 
Electrode  Tip  Radius  and  Spacing  from 
Aircraft. 
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explanation  for  the  equality  in  view  of  the  difference  in 
y  and  yS*  fp  is  that  }f  is  decreased  for  the  negative  elec¬ 
trode  by  the  effects  of  surface  absorption  of  oxygen  be¬ 
coming  fortuitously  equal  to^fp  for  air. 

A  difference  in  the  magnitudes  of  positive  and 
negative  corona  current  for  the  same  electrode  voltage 
has  been  observed  (Reference  5).  This  is  caused  by  a 
difference  in  the  clearing  time  of  the  choke-off  space 
charge  which  affects  the  pulse  shape  and,  therefore,  the 
average  current.  The  rate  of  clearing  is  approximately 
twice  as  rapid  with  the  negative  electrode  since  the 
electrons  traverse  part  of  the  gap  as  electrons  and  attach 
to  ions  later.  Also,  the  mobility  of  negative  ions  in 
air  is  30  per  cent  greater  than  that  of  positive  ions 
due  to  its  smaller  mass  to  area  ratio.  The  number  of 
ions,  m,  in  streamers  and  burst  pulses  is  of  the  same  order 
of  magnitude  as  in  the  negative  Trichel  pulses.  The 
pulse  repetition  rate,  Fr,  however,  is  reduced  some  50  per¬ 
cent  for  positive  corona  current  because  of  the  longer 
clearing  time.  Consequently,  near  onset  for  the  same 
conditions,  the  negative  corona  currents  in  air  are  of  the 
order  of  twice  as  great  as  the  positive  corona  currents. 

For  both  polarities  an  Ohm’s  law  regime  exists  near  onset. 

CORONA  DISCHARGING  CURRENT  IN  STILL  AIR 

An  equation  for  the  magnitude  of  the  corona  dis¬ 
charging  current  is  presented  in  Reference  2,  page  32,  as 
will  be  further  discussed  in  Section  D  (equation  D-50). 

In  still  air,  this  equation  reduces  to 

X  ■  ACV*  (^-  '/c)  (C-22) 

where  A  and  C  are  constants  for  a  particular  ccrona  probe- 
fuselage  configuration.  Figure  5  presents  experimental  data 
of  corona  current  at  zero  airspeed  measured  during  Contract 
DA  44-177-AMC-T(3)  and  from  other  data  obtained  by  this 
contractor.  For  a  value  of  AC  =  0.83  x  lO"^  and  Vc  =  2000 
volts,  equation  (C-22)  provides  the  following  comparison 
with  the  experimental  data. 
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A  Data  from  Reference  21 
A  Unpublished  Data 


Applied  Voltage,  Va  x  10^ 


Figure  5.  Corona  Discharging  Current  in  Still  Air. 
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TABLE  I 


Theoretical 


Experimental 


Va  x  10“ 3  Volts 

I  x  10"6 

amps 

I  x  10’ ^  amps 

0 

0 

0 

1 

-0.0083 

(maximum 

negative 

value) 

0 

2 

0 

0 

10 

0.664 

1 

23 

4.0 

4.0 

100 

81.5 

78.0 

It  should  be  noted  that  whereas  experimentally,  and 
obviously  physically,  the  current  is  zero  up  to  the  corona 
threshold  Vc  =  2000  volts,  the  theoretical  expression  shows 
negative  currents  for  values  of  Va  <  Vc. 

The  reason  for  this  discrepancy  is  that  the  assump¬ 
tions  made  in  formulating  equation  (C- 22)  do  not  consider 
the  discontinuity  that  the  threshold  voltage  represents  in 
the  basic  physical  process  of  corona  action.  The  discrepancy, 
however,  is  seen  to  result  in  a  practically  insignificant 
small  value  of  negative  current. 
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SECTION  D:  CALCULATION  OF  CORONA  CURRENT  IN  MOVING  AIR 


Derivation  of  an  expression  for  corona  current  serves 
to  make  evident  the  various  factors  influencing  discharge 
current  magnitude.  This  information  enables  design  of  the 
corona  electrode  for  maximum  discharge  current  capability. 

Such  optimization  makes  possible  light,  small  discharger 
systems.  The  derivation  is  based  on  the  corona  theory 
given  in  the  preceeding  section.  This  section  first  gives 
a  general  description  of  the  current  discharge  from  a  corona 
electrode  mounted  on  a  helicopter  under  in-flight  conditions. 
Detailed  descriptions  of  the  factors  influencing  the  discharge 
follows.  The  section  includes  a  derivation  of  an  equation 
for  corona  current  involving  these  factors.  The  ion  tra¬ 
jectory  is  also  stated. 

AIRCRAFT  CHARGE  DISSIPATION  FROM  THE  CORONA  El ECTRODE 

The  discharger  system  applies  a  large  voltage 
difference  between  the  corona  electrode  and  the  aircraft 
surface  in  order  to  create  large  electrical  field  strength 
at  the  electrode.  This  field  strength  causes  ionization 
of  the  surrounding  atmosphere  and,  consequently,  causes 
the  presence  of  positive  ions  and  electrons  near  the 
electrode  (Section  C) .  Since  the  polarity  of  the  elec¬ 
trode  is  that  of  the  electrostatic  charge  picked  up  by 
the  aircraft,  the  particles  (positive  ions  or  electrons) 
attracted  to  the  electrode  neutralize  the  aircraft  charge. 

The  trajectory  of  the  repulsed  particles  is  of  importance 
for  the  following  reasons.  Any  accumulation  of  repulsed 
particles  near  the  corona  electrode  constitutes  a  space 
charge  which  limits  the  magnitude  of  discharge  current 
and  reduces  the  ionizing  electric  field.  In  addition, 
any  return  of  repulsed  particles  to  the  aircraft  reduces 
discharger  system  effectiveness  by  reducing  the  net  dis¬ 
charge  current.  The  trajectory  of  the  repulsed  particles 
is  determined  by  the  total  force  exerted  on  them.  Because 
of  their  small  mass,  the  ions  will  travel  in  the  direction 
of  the  total  vector  force  at  any  point.  This  force  is  the 
vector  sum  of  the  electrical  force  and  the  viscous  air- 
stream  or  wind  force  on  the  particles  at  any  instant.  That 
is, 
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0>-l) 


where  Ft 


vector  total  force  on  the  repulsed 
particles  at  any  instant, 


electric  fores  vector  on  the  repulsed 
particles,  and 


Fw  =  vector  viscous  wind  force  at  any  instant. 


The  two  components  of  the  total  force  are  discussed 
in  detail  in  this  section. 

In  summary,  the  phenomena  determining  the  magnitude 
of  corona  discharge  current  and  which,  therefore,  must  be 
considered  in  the  derivation  are  as  follows: 

1.  Electric  field  of  specific  electrode 
geometry  influenced  by  a  distribution 
of  space  charge. 

2.  Ionization  of  air  by  this  electrical 
field. 

3.  Space  charge  density  buildup  resulting 
from  this  ionization  as  a  function  of 
time  and  position. 

4.  Effect  of  wind  and  electric  field  on 
space  distribution  and  ion  trajectories. 

EFFECT  OF  AIRCRAFT  ON  DISCHARGER  ELECTRICAL  FIELD 

The  e*. metrical  component,  in  equation  (D-l) 
is  significantly  affected  by  the  electrically  conductive 
objects  in  proximity  to  the  corona  electrode.  The  dis¬ 
charger  system  can  be  considered  to  operate  completely 
isolated  from  all  objects  except  the  aircraft  structure 
since  the  discharger  functions  only  during  flight.  The 
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proximity  of  a  conductive  object  connected  to  the  airframe 
will  increase  the  recirculated  current  percentage,  and  hence 
decrease  the  probe  effectiveness. 

The  field  distribution  near  a  corona  probe  has  been 
analyzed  in  Section  B.  The  analysis  did  not  include  the 
effect  of  space  charge  density  due  to  the  motion  of  ionized 
particles.  This  density  is,  however,  very  small  except  in 
the  zone  very  near  to  the  corona  probe,  and  in  this  zone, 
the  effect  of  the  space  charge  was  discussed  in  Section  C. 
Hence,  the  model  of  the  corona  discharge  in  air  is  a  blend 
between  the  corona  discharging  phenomenon  taking  place  in 
the  zone  near  the  probe  followed  by  an  ionic  distribution 
determined  by  the  electrostatic  and  viscous  forces  pre¬ 
dominant  in  zones  relatively  far  away  from  the  probe  tip. 

The  current  distribution  in  these  zones  will  be  analyzed 
later  in  this  section. 

COMPONENT  OF  ION  MOTION  DUE  TO  ELECTRICAL  FIEUD  OF  DISCHARGER 

It  is  instructive  to  compare  ion  motion  under 
actual  discharger  operating  conditions,  that  is,  in  the 
atmosphere,  with  their  motion  in  a  vacuum.  In  the  latter 
case  the  ions  travel  along  the  electric  field  lines  with¬ 
out  interference  from  the  air  molecules  with  an  increasing 
velocity  v  given  by  Newton’s  second  law.  The  differential 
equation  which  must  be  solved  to  obtain  v  is 

Ft  *  (D-2) 

is  force  vector  on  ion  equation  (B-9), 
is  ion  mass, 

is  vector  ion  velocity,  and 

is  vector  electric  field  caused  by  dis¬ 
charger  system,  as  shown  in  equation  (B-7). 

The  presence  of  the  atmosphere,  however,  slows 
up  the  motion  of  ions  in  an  electric  field  because  of  the 
repeated  collisions  with  neutral  molecules  which  result 
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where 


re 

mi 

v* 

t 


in  kinetic  energy  losses.  At  each  impact  the  ion  loses 
much  of  its  velocity  in  the  field  direction.  After  a 
large  number  of  collisions  have  occurred,  the  ion  has 
an  average  or  drift  velocity  in  the  field  direction  which 
is  directly  proportional  to  and  inversely  proportional 
to  the  atmospheric  density.  This  drift  velocity  is 
constant  for  a  given  E  rather  than  increasing  as  it  does 
in  vacuum.  The  drifc  velocity  divided  by  the  electric 
field  is  therefore  a  constant.  This  ratio  is  termed 
"mobility"  and  depends  only  on  gaseous  composition  and 
density  (Reference  6).  Thus,  in  still  air,  ions  in  the 
atmosphere  still  travel  along  the  lines  of  force  but  with 
a  constant  velocity.  Assuming  zero  initial  velocity, 
the  expression  for  ion  velocity  in  atmosphere  is 


V  * 


(D-3) 


where  v* 

ion  velocity, 

ki  = 

mobility  constant  of  ions,  and 

E 

field  strength. 

Since  the  electric  field  lines  of  the  discharger  system 
extend  between  the  corona  electrode  and  the  aircraft 
surface,  there  is  a  definite  tendency  in  still  air  for 
discharge  current  to  return  to  the  aircraft.  If  it  were 
not  for  the  airstream  past  the  discharger  electrode, 
special  measures,  such  as  the  application  of  insulative 
coatings  to  the  aircraft,  would  have  to  be  considered. 
The  wind  prevents  return  current  by  exerting  a  force  on 
the  ions  which  diverts  them  from  the  electric  field 
lines. 

COMPONENT  OF  ION  MOTION  DUE  TO  AIRSTREAM  PAST  DISCHARGER 
ELECTRODE 


The  ion  motion  caused  by  the  airstream  is  com¬ 
pletely  determined  if  the  viscous  wind  force  on  the  ion 
and  the  initial  ion  velocity  are  known.  If  the  viscous 
force  on  the  ions  is  much  larger  than  the  electrical 
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force  their  trajectory  is  along  the  lines  of  viscous 
force;  if  not,  the  ion  trajectory  is  along  the  vector 
resultant  of  the  viscous  and  electrical  force. 


The  magnitude  of  viscous  force  on  a  body  in  an 
airstream  is  given  by  conventional  aerodynamic  theory  as 
a  function  of  the  coefficient  of  drag.  This  coefficient 
is  a  function  of  the  shape  of  the  object  and  the  reynolds 
number  of  the  airstream.  Reference  10  suggests  this 
technique  can  be  used  to  determine  the  magnitude  of 
viscous  force  on  the  ions  of  the  corona  discharge  current 
if  the  ion  is  considered  a  small  rigid  sphere.  This 
reference  gives  drag  coefficient  information  for  spheres 
with  diameter  as  small  as  2  x  10“^  centimeter.  The  dia¬ 
meter  of  the  ions  of  the  discharge  current  are  much 
smaller  than  this  value  when  first  formed.  However,  rapid 
attachment  of  these  ions  by  larger  particles  of  the  air 
augments  the  particle  diameter  by  several  orders  of 
magnitude.  Thus,  the  use  of  formulas  in  Reference  10  to 
calculate  an  approximate  drag  coefficient  for  the  ions  of 
the  discharge  current  is  a  justified  extrapolation. 


Accordingly,  the  viscous  force  exerted  by  the 
airstream  on  the  ions,  f  the  aircraft  discharging  current 
is  given  by 

Fw  *  Op  <%h  S  (D-4) 


where  qr  =  dynamic  pressure  caused  by  airst  *eam, 

S  =  ion  projected  area  (assumed  spherical), 

Op  =  coefficient  of  drag;  and 

Fw  =  viscous  airstream  force. 

For  practical  use  of  equation  (D-4),  equations  for  the 
quantities  appearing  in  it  must  be  given.  These  are: 


lr  =  1/2  (A  *x) 


(D-5) 


CD  =  27/  (Re)'84 


(D-6) 
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Rt  * 


(D-7) 


wd 

Y 

Applying  these  relationships  to  equation  (D-4)  gives  a 
more  useful  expression  for  viscous  force.  This  is 

K/m  Vx  Ce  /> A  WXS*  339/TV'*¥(WD)  fiA  0>“8) 

In  the  preceeding  four  equations. 

/>A  is  airstream  density  (free  stream) , 

W  is  airstream  speed, 

Re  is  reynolds  number, 

d  is  diameter  of  ion,  and 

V  is  kinematic  viscosity  of  airstream. 

Equation  (D-8)  displays  the  dependence  of  vis¬ 
cous  force  7^  on  environmental  conditions.  It  is  seen 
that  viscous  force  increases  nearly  linearly  with  increase 
of  airstream  speed  and  kinematic  viscosity.  The  eifects 
of  temperature  and  pressure  on  viscous  force  nre  contained 
in  the  direct  proportionality  of  to  free  stream  density 
/’A-  Since  /»A  ls  directly  proportional  to  ambient 
pressure  and  inversely  proportional  to  temperature,  vis¬ 
cous  force,  7^,  increases  as  ambient  pressure  increases  if 
temperature  remains  constant.  Also  the  viscous  force 
decreases  with  increase  cf  atmospheric  temperature  if 
ambient  pressure  remains  constant. 

ION  TRAJECTORY  DUE  TO  VISCOUS  AND  ELECTRICAL  FORCE  OF  DIS¬ 
CHARGER  SYSTEM 

By  use  of  equations  (D-l),  (B-9),  and  (D-8),  the 
forces  acting  on  an  ion  can  be  computed. 
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0.f¥ 


(wo) 


(D-9) 


typical  values  obtained  from  equation  (D-9)  for 
Fe  and  Fw  are  shown  in  Figure  6  and  Figure  7  for  h=c  and 
h  =  c/2,  respectively.  In  these  figures,  the  following 
values  were  utilized: 


Rt  =  0.002  inches 

Vp  =  50,000  volts 

_  Q 

q  =  1.6  x  10  coulombs 

c  =24  inches 
d  =  1.43  x  10“8  inches 

W  =  720  inches/second 

As  seen  from  Figures  6  and  7,  near  the  corona  elec¬ 
trode  the  electric  component  of  force,  Fe,  on  the  ions  is 
much  larger  in  magnitude  than  the  viscous  airstream  com¬ 
ponent  Fw.  Thus,  the  ion  trajectory  near  the  electrode 
will  be  in  the  direction  of  the  electric^f ield  lines,  that 
is,  the  trajectory  has  the  direction  of  Fe  at  points  near 
the  corona  electrode. 
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Electric  and  Viscous  Forces  Acting  on  Charged  Particle 


r , 


Figure  6.  Force  Vectors  Acting  on  a  Charged  Particle 
at  h  =  c. 
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Electric  and  Viscous  Forces  Acting  on  Charged  Partic 


Vertical  Distance  Below  Corona  Point,  R,  Inches 


Figure  7.  Force  Vectors  Acting  on  a  Charged  Particle 
at  h  =  c/2. 
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of  the  fetal  %rce*  5$,  becomes  quite  small  in  coepaxM^sa 
m  the  vistas  airs  treats  component  which  remains  relatively 
constant  hence,  ^  has  essentially  the  direction  ol 

the  doww^^  Velocity  factor, 

EFFECT  OF  ElfemODE  TIP;  RADIUS 

Fo$  a  discharge  electrode  of  point  geometry,  there 
is  an  optimum  radius  of  curvature  which  represents  a  com¬ 
promise  between  the  factors  described  below.  Because  of 
deterioration,  however,  it  is  impossible  to  maintain  this 
optimum  radius .  The  operating  conditions  determined  by  tip 
radius  and  tfcf  significance  of  these  conditions  to  system 
perfprmancf  are  discussed  in  the  following  paragraphs. 

Magnitude  or  Field  Intensity  in  the  Vicinity  of  the/  Corona 
Electrode  it  a  Given  Voltage 

The  field  strength  for  given  applied  voltage  and 
ambient  charge  density  about  the  corona  electrode  is  com¬ 
pletely  determined  by  the  electrode  configuration.  Corona 
current  oeftijra  only  when. the -field  strength  is  equal  to  or 
greater  tWi  a  value  vi  approxiisa-wly  uU  ktiovolto  >.  er 
centimeter  (Reference  19) .  The  cprona,  current  pulse  re¬ 
petition  rate,  Fr,  Increases  as  the  field  strength  in¬ 
creases  and,  therefore,  according  to  equation  (C-17),  the 
corona  current  also  increases  with  field  strength*  Small 
radius  of  curvature  of  the  corona  electrode  tip  will 
minimize  the  voltage  required  to  produce  a  given  field 
strength  since  the  field  increases  as  the  radius  of 
curvature  of  the  tip  decreases.  The  voltage  required  for 
corona  onset  depends  on  the  electrode  tip  radius  in  the 
manner  shown  in  Figures  3  and  4.  Magnitude  of  corona 
current  increases  as  the  corona  onset  potential  decreases 
in  accordance  with  equation  (D-50). 

Extent  of  the  High  Field  Region  Around  the  Electrode  Tip 

As  shown  in  equation  (C-17),  the  corona  current 
level  is  proportional  to  the  number  of  ions  in  one  pulse. 

A  large  value  of  m  will  result  only  if  the  extent  of  high 
field  region  about  the  corona  electrode  is  appreciable. 

The  extent  of  high  field  region  increases  as  the  radius 
of  curvature  of  the  tip  decreases. 
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Since  current  density  in  corona  discharge  is  constant 
over  a  vide  range  of  voltage,  the  current  is  dir)/1  tly  pro¬ 
portional  to  the  point  surface  area,  too  small  a  point  is 
very  prone  to  performance  degradation  due  t;o  positive  ion 
bombardment  (sputtering) .  Electrode  tip  deterioration  by 
sputtering  results  in  a  tip  with  a  much  increased  radius 
of  curvature  and,  consequently,  an  increased  surface  area, 
(Reference  18) .  Too  large  a  surface  area  may  mean  a  loss 
in  field  intensity. 

Effect  of  Wind 


The  presence  of  wind  produces  viscous  forces  on 
the  ions.  In  the  close  vicinity  of  the  probe  tip,  these 
forces  are  substantially  smaller  in  magnitude  than  the 
electrostatic  forces,  which  in  that  zone  are  necessarily 
very  high  to  produce  ionization.  Hence,  the  resultant 
force  vector  is  essentially  equal  to  the  electrostatic 
force  vector.  This  implies  that  the  effects  on  the 
probe  performance  of  tip  radius  and  the  airspeed  are 
essentially  negligible  Vector  .ret-.' o  of  tM 

forces  acting  on-  the  Idris  in  the  vicinity  of  the  corona 
tip  is  shown  in  the  sketch  below. 
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The  theoretical  procedure  for  formulatiers  of  the 
system  of  equations  which  would  he&e  to  be  solved  for 
determination  of  the  exact  expression  for  corona  curfefit  is 
now  discussed.  The  resulting  equations  are  very  complex 
and  exact  mathematical  solutions  require  techniques  which 
are  beyond  the  scope  of  this  program.  Formulation  of  the 
equations,  however,  does  have  the  value  of  clarifying  the 
physical  mechanisms  which  occur  and  thus  provide  guidelines 
for  sound  corona  prpbe  design  and  experimental  optimizing 
modifications . 

The  formulation  of  the  equations  is  achieved  by 
mathematically  adding  up  the  changes  in  the  number  of 
charged  particles  in  an  infinitesimal  volume  at  a  given 
location  in  the  corona  discharge  during  an  infinitesimal 
time  interval  due  to  the  various  physical  mechanisms  to 
obtain  the  overall  rate  of  change  for  each  particle.  Since 
there  are  three  different  particles  involved  (electrons, 
posltiya  ions  ,  and '  negative  ions  due,  to  electron  attachment 
to  oxygen  molecules) ,  three  equations  are  necessary.  The 
fourth  and  final  equation  is  the  well-known  poisson  electro¬ 
static  equation  relating  the  number  of  charged  particles  per 
unit  volume  to  electrical  field.  An  exact  solution  for  this 
system  of  equations  would  give  the  concentration  of  each  of 
the  particles  and  the  electrical  field  as  a  function  of  time 
and  location. 

The  physical  mechanisms  described  in  the  equations, 
a  particle  population  balance  diagram  for  an  infinitesimal 
volume,  and  derivation  of  the  equations  are  given  below. 

A  general  expression  showing  how  the  value  of  corona  current 
could,  in  theory,  be  obtained  from  the  solution  of  the  above 
four  equations  will  also  be  given. 

Assumptions  are  made  as  described  in  the  following 
derivation. 
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1.  For  egctffo  4&84ttoa  generated  by-  Ionization, 
a  .pcf&iv&xte  loti  Is  generated  , 

2.  Positive  oxygen  and  nitrogen  i^ns  have 
the  &diwe  diffusion,  recombination,  and 
mobility  constants . 

3.  Posit ions  do  not  produce  ionization 
Because  of  their  relatively  large  mass. 

4.  Although  electron  and  positive  ions  are 
generated  by  the  same  ionlzatinp  event, 

:lsi  not  equal  n£  because  of  ?:i .jctron 
Qttacilmpnt  to  form  negative  icr_.', . 

5.  No  combination  of  positive  and  negative 
ions  occurs. 

6„  Molecules  of  the  atmosphere  ar<*  fir^le 
lionized.  w  -  - 

The  physical  mechanisms  of  importance  in  determining 
electrostatic  discharger  system  corona  current  follow. 

Particle  Generation  Due  to  Ionization  By  Electrons 

In  advancing  a  unit  distance  in  the  atmosphere  in 
the  direction  of  the  electrical  field,  a  single  electron 
creates  «*<  (Townsend’s  first  ionization  coefficient  dis¬ 
cussed  in  Section  C,  page  13  )  new  electrons  and  positive 
ions.  The  Increase  in  charged  particles  dn  caused  by  ni 
electrons  in  advancing  a  distance  dx,  is  given  by  the  re¬ 
lation 


<//!  »  ci.  f) 

(D-10) 

The  rate  of  electron  increase  due  to  this  mechanism  in  the 
infinitesimal  volume  (of  unit  cross  section)  is 


45 


nt  v 
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1ft  * 


By  using  equation  (D-3) ,  there  follows 


oCrt, 

ct-t 


«  oi.  n,k ,  | e, 


(D-12) 


where  m  is  number  Of  electrons  located  within  the  in 
finitesimal  volume,  that  is,  electron  concentration. 


If  the  assumption  is  made  that  one  positive  ion  re 
suits  from  each  electron  generated  in  this  mechanism,  the 
contribution  to  the  positive  ion  population  due  to  ioni¬ 
zation  by  electrons  is 


At 

d-t 


CD-13) 


where  n^  is  the  number  of  positive  ions  in  the  infinitesi¬ 
mal  volume.  The  viscous  forces  acting  on  the  electron  can 
be  considered  in  the  preceding  two  equations  by  use  of  an 
’’effective”  electric  field.  This  is  accomplished  by  assum¬ 
ing  E  to  consist  of  the  actual  electric  field,  Ej,  plus  an 
equivalent  electric  field*  which  expresses  tne  viscous 
effects  on,  the  particular  ion.  That  is 


E  =  tj,  +  \  (D-14) 

where  E^  =  T^/e 

^w  =  9[)qrS/e 

Since  S  is  very  small  for  electrons,  the  effects  of 
viscous  forces  in  contributing  to  ionization  are  negligible 
in  comparison  to  the  effects  of  electric  forces.  Viscous 
forces  are  important  in  determining  distribution  of  ions 
through  other  mechanisms  as  described  below.  Ionization 
by  electrons  is  not  a  factor  in  causing  negative  ions. 
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Particle  Generation  Due  to  Volume  Ionization 

Here,  reference  is  Trade  to  the  contribution  made  by 
naturally  occurring  ambient  ionization  normally  always  pre¬ 
sent  and  due  to  cosmic  rays  and  natural  radioactivity. 

This  mechanism  is  represented  mathematically  by  the  symbol 
qi,  where  equals  natural  time  rate  of  ionization.  The 
term  q^  is  the  rate  of  positive  ion  generation  as  well  as 
electron  generation  through  natural  causes.  This  stems 
from  the  fact  that  for  each  electron  produced,  one  positive 
ion  is  produced,  that  is,  the  molecules  of  the  atmosphere 
are  being  single  ionized  by  these  natural  causes.  Volume 
ionization  does  not  contribute  to  negative  ion  generation. 

Particle  Depletion  Due  to  Recombination 

In  this  mechanism,  electrons  reunite  with  positive 
ions  causing  the  loss  of  the  two  particles  from  their  re¬ 
spective  populations  and  the  possible  emission  of  light  in 
the  process.  The  number  of  individual  recombinations  during 
an  incremental  time,  dt,  for  a  region  in  which  n^  and  ^ 
are  uniform  concentrations  of  particles  per  unit  volume  is 

dn  •  d*  (d-15) 

where,  by  definition,  *<r,  is  the  coefficient  of  recom¬ 
bination  (Reference  19,  page  20).  The  quantity  *<r  de¬ 
pends  on  pressure  and  temperature  since  these  conditions 
affect  the  number  of  positive  and  negative  ions  present 
and  available  for  recombination  in  a  given  volume  of  air. 

For  air  at  0°C  and  760  ram  pressure,  «<r  is  1.71  x  10'® 
cubic  centimeter/ second  (Reference  6,  page  99). 

Particle  Depletion  Due  to  Diffusion 

Wherever  the  concentration  of  charged  particles 
varies  from  point  to  point  in  a  gas,  that  is,  where  a  con¬ 
centration  gradient  exists,  there  will  be  a  flow  of  these 
particles  from  regions  of  high  concentration  to  regions  of 
lower  concentration.  Thus,  diffusion  produces  a  de-ionizing 
effect  in  high  concentrated  regions  and  an  ionizing  effect 
in  the  adjacent  regions  of  lower  ion  concentration.  If 
n^  is  the  number  of  charged  particles  at  a  given  point,  the 
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concentration  gradient  is  S?n*  where  S7  is  the  differential 
operator  of  vector  calculus,  the  rate  of  increase  or  de¬ 
crease  into  or  out  of  the  small  volume  being  considered  in 
the  particle  population  balance  is  directly  proportional  to 
the  divergence  of  the  concentration  gradient,  that 
or  Thus  rate  of  particle  increase  or  decrease  due 

to  diffusion  is  given  mathematically  by 


Uf)c 

lux 


Dc  V2/it‘ 


(for  i  =  1,  2,  and  3) 


(D-16) 


where  D£  is  the  constant  of  proportionality  or  diffusion 
constant. 


The  value  of  dn/dt  may  be  positive  or  negative  corresponding 
to  ionization  or  de-ionization  (depletion),  respectively, 
according  to  whether  the  net  gradient  over  the  surfaces  of 
the  infinitesimal  volume  (the  vector  divergence)  is  positive 
or  negative  (Reference  11,  page  189).  The  diffusion  coef¬ 
ficient  is  related  to  the  mobility  constant  for  a 
given  ion  by 


Oi  * 


kcp 

*re 


(D-17) 


where  N  is  the  number  of  molecules  per  cubic  centimeter 
at  pressure  p  (Reference  19,  page  30). 


The  ion  mobility  constant  is  discussed  below  in 
connection  with  the  physical  mechanism  of  particle  gener¬ 
ation  or  depletion  due  to  electrical  charge  motion  in¬ 
duced  by  the  electric  field.  is  inversely  proportional 
to  pressure  and,  for  atmosphere  with  water  content,  directly 
proportional  to  the  1.75  power  of  temperature. 


Negative  Ion  Generation  Due  to  Electron  Attachment 

As  mentioned  in  the  description  of  the  events 
occurring  in  corona  discharge  (Section  B),  the  attachment 
phenomena  causes  the  pulsed  nature  of  the  discharge.  In 
this  phenomenon,  electrons  are  attached  by  neutral  molecules 
of  the  oxygen  of  the  atmosphere  when  these  electrons  reach 
regions  of  the  discharge  where  E/p  has  been  reduced  to  20 
volts /centimeter/  nin  Hg  or  less  (Reference  19,  page  28). 

The  coefficient  of  attachment,  <f  ,  which  is  the  average 
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number  of  collisions  the  electron  srust  make  before  attaching 
itself  to  a  neutral  atom,  is  not  related  to  any  property  of 
the  atmosphere,  except  to  its  electronegative  (tendency  to 
capture  electrons  due  to  0£  content)  character.  If  p  is 
the  mean  time  for  an  electron  to  become  attached  to  a  mole¬ 
cule  of  the  atmosphere,  then  the  rate  of  electron  depletion 
from  the  infinitesimal  volume,  or  alternatively,  the  rate 
of  negative  ion  generation,  is  given  by 

—p—  -  JU  (D-18) 

d*  Y 

The  value  of  Y  for  air  at  0°C  and  760  mm  Hg  is  6.3  x  10”  ^ 
seconds  (Reference  6,  page  97).  There  is  a  dependence  of 
Y  on  temperature  and  pressure,  since  p  equals  the  coef¬ 
ficient  of  attachment  divided  by  the  number  of  electron 
collisions  per  second  with  air  molecules,  the  latter  being 
temperature  and  pressure  dependent. 


Particle  Generation  or  Depletion  Due  to  Electrical  Charge 
Motion  Induced  by  the  Electric  Field 


From  basic  electric  theory  (Reference  8,  page  187) 
the  rate  of  charge  buildup  or  decrease  (depending  on  alge¬ 
braic  sign  of  the  final  result),  dq/dt,  in  the  infinitesimal 
volume  under  consideration  is  given 


ff  J;  •  aCs** (for  i  =  1,  2,  3)  (D-19) 

where  is  current  density  of  charged  particles,  and  d!f 
Is  infinitesimal  area  element.  a 

If  represents  the  volume  charge  density  at  the 
location  of  the  Infinitesimal  volume  under  consideration 

' ^7 ' '  7*ffL  A  ff  (D‘20) 


Usinfc  stokes  theorem, 


![sfc * II[*'  f‘  s*  <°-2i> 

where  the  negative  sign  and  differentiation  with  respect  to 
time  has  been  brought  under  the  triple  integral  sign.  Thus, 
the  following  equation  for  rate  of  change  of  volume  charge 
density  is  obtained  by  equating  integrands  in  the  preceding 
equal  triple  integrals: 


V  *  (D-22) 

The  preceding  equation  will  be  applied  to  each  of 

the  three  types  of  charged  particle  involved  in  the  corona 
discharge  in  deriving  the  discharge  equations.  For  electrons, 

:  -Y),eV,  (D-23) 

The  minus  sign  in  the  preceding  equation  is  due  to  conven¬ 
tion.  Also, 

J°t  *  h,e  (D-24) 

Thus,  the  rate  of  change  of  electron  population  in  the 
infinitesimal  volume  due  to  field  induced  charge  motion, 
dn^/dt,  is  given  by 

-e  v,  -  -  e  (D-25) 

at 

or 
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The  electron  velocity  can  be  replaced  immediately 
by  one  of  the  quantities  whictef^he  system  of  equations  to 
be  derived  should  involve  and,  also,  give  knowledge  of. 
This  quantity  is  the  electric . field  intensity  E.  The 
relationship  to  be  used  to  enable  this  replacement  is 
equation  (D- 3).  Thus, 

-  v>  n,  M,  /T  (d-26) 

Again  can  be  considered  to  consist  as  the  sum  of  the 
actual  electric  field  and  an  equivalent  electric  field  to 
represent  the  viscous  force  effects. 

The  mobility  constant  for  ions,  k^,  is  independent 
of  E/p  over  a  wide  range  of  E/p  values  provided  that  the 
speed  of  drift  in  the  field  is  appreciably  less  than  the 
velocity  of  thermal  agitation.  The  product  of  .electron 
mobility  and  pressure,  k^  x  p,  in  air  is  approximately 
constant  from  an  E/p  value  of  0  to  22  volts /centimeter/mm 

Hg.  The  approximate  value  of  the  ki  x  p  product  is 

cm 

5  x  10"^  sec  111111  (Reference  19,  page  35.) 

V/cm 

For  positive  ions, 

*  n2eV2  (D-27) 


where  in  fehis  case,  V2  is  the  positive  ion  velocity  and  e 
is  positive  ion  charge  (equal  to  that  of  an  electron  because 
the  air  molecule  has  been  deprived  of  one  electron  in  form¬ 
ing  the  positive  ion) . 

The  signs  in  the  preceding  equation  are  such  that 
positive  ion  current  density  is  in  the  same  direction 
(not  the  opposite  direction  as  in  the  case  of  electrons) 
as  the  ion  velocity  in  accordance  with  convention.  Also, 


A  *  n*e 


(D-28) 
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Thus,  the  rate  of  change  of  positive  ion  population  in  the 
infinitesimal  volume  due  to  field  induced  charge  motion, 
d^/dt,  is  given  by 

c//ix 

V**zv*~e  —  (D-29) 


oCt 


or 


/>XVZ 


*  - 


ctnx 

oC< 


The  positive  ion  velocity  can  be  replaced  as  in  the  case  of 
electrons  by  the  total  electric  field  intensity  eL  experienced 
by  the  positive  ion.  contains  an  electric  field  component 
equivalent  to  the  air  stream  viscous  forces  on  positive  ions. 
For  positive  ions,  this  viscous  force  term  may  not  be  negli¬ 
gible.  Using  the  ion  mobility  constant  for  positive  ions, 


Thus, 


vx  *  kx  <D"30> 

~UTC  **  (D-31) 


The  mobility  constant  for  atmospheric  (dry)  posi¬ 
tive  ions  (singly  charged)  at  0°C  and  760  mm  Hg  is  1.36 
centimeters /second/volt /centimeters  (Reference  6,  page  38). 
For  negative  ions  (those  resulting  from  the  attachment  of 
electrons  by  air  molecules  due  to  their  electronegative 
property  and  which,  by  causing  the  pulsed  nature  of  corona 
discharge  as  described  elsewhere  in  this  report,  play  a 
significant  role  in  electrostatic  discharger  operation) , 


(D-32) 

Again,  the  negative  sign  in  the  preceding  equation  iV  due 
to  the  definition  of  current  density  in  terms  of  positive 
charge  where  here  negative  charge  is  being  considered.  Also, 


A  ■ 


(D-33) 
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Thus,  the  rate  of  change  of  electron  population  in  the  in¬ 
finitesimal  volume  due  to  field  induced  charge  motion,  d^/dt 
is  given  by 

-*  -e  (D-34) 


or 


■  #■ 


Replacing  V3  by  k3E*3,  as  before^  where  IC3  is  the  mobility 
constant  for  negative  ions  and  E3  is  the  total  electric  field 
(contains  a  component  6^3,  equaling  the  electrical  equiva¬ 
lent  of  the  airstream  viscous  forces  on  the  negative  ion) 
gives 


dn9 


=  V*  K* 


(D-35) 


For  negative  ions  the  viscous  force  term  may  not  be 
negligible.  The  mobility  constant  for  negative  ions  of  in¬ 
termediate  size  (smaller  than  "Langevin"  ions)  at  0°C  and 
760  mm  Hg  is  of  the  order  of  7  x  10’ ^  square  centimeter/ 
volt  second  (Reference  11,  page  172). 


For  still  smaller  singly  charged  negative  ions 
(nonclustered)  at  0°C  and  760  mm  Hg,  the  mobility  constant 
is  2.1  square  centimeter/volt  second  (Reference  6,  page  38).. 
Figure  8  shows  a  gaseous  ion  particle  balance  for  the  in¬ 
finitesimal  volume  under  consideration  applying  all  the 
ionization  and  de-ionization  mechanisms  described  above. 
With  the  information  summed  up  in  this  diagram,  three 
equations  involving  the  concentrations  ni ,  ti^  of  the 
three  ions  in  the  corona  discharge  will  be  written  below. 

A  fourth  equation  must  be  obtained  which,  in  conjunction 
with  the  aforementioned  three,  will  enable  solution  in 
theory  for  the  four  quantities  n^,  n£,  n^,  Ep  at  every 
point  in  the  discharge.  This  equation  is  basically  the 
well  known  poisson  equation  and  can  be  derived  from  the 
gauss  law  of  electrostatics  (Reference  8,  page  43)  and 
vector  analysis.  This  law  states  that  the  normal  com¬ 
ponent  of  the  electric  displacement  vector,  taken 
over  the  surface  of  any  closed  volume  equals  the  net 
charge  enclosed  in  that  volume.  As  an  equation,  the 
gauss  law  takes  the  form 
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Figure  8*  Particle  Balance  Diagram  for 
Infinitesimal  Volume. 
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where 


K  *  ‘ds*'2  ‘//La 

**a 


(D-36) 

(for  1=1,  2,  and  3) 


dS.  is  an  element  of  surface  area, 

q  is  net  charge  enclosed  in  the  infinitesimal 
volume , 


/>  i  is  net  charge  volume  density  of  particles,  and 
dVa  is  an  element  of  volume. 


Using  the  divergence  theorem,  the  preceding  equation  leads 
to 

ff  Pc  ■  At  *  [[[' V-  f  f  ffi  (D-37) 

'r*  %.  iL 

Since  the  limits  of  integration  in  the  second  and 
third  members  of  the  preceding  equation  are  the  same,  the 
integrands  must  be  equal.  That  is, 

(D-38) 

Electric  theory  gives  a  relationship  by  which  f  may 
be  replaced  in  the  preceding  equation  with  6ne  of  the  quanti¬ 
ties  to  be  found  with  the  four  equations,  that  is,  Ep.  This 
relationship  is 

P<  *  €.  tr  (d-39) 

>$'■  ' 

where  €0  is  electric  permittivity  of  free  space  having 
a  value  8.85  x  10‘12  farad/meter  in  the  MKS  system. 


In  addition, the  ion  concentrations  (number  of  ions 
per  unit  volume  at  an  arbitrary  location  of  the  infinitesi¬ 
mal  volume)  nj,  n2>  n^  for  the  three  ions  of  the  discharge 
for  which  the  equations  are  to  be  solved  can  be  put  in  the 
equation  by  replacement  of  the  quantity  JD  .  The  relation¬ 
ship  through  which  this  is  accomplished  is 


f  *  A  *  /’a  *  J°,  *  ('>,+'>*)] 


(D-40) 
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where  e  is  the  electronic  charge  with  a  value  pf  1.6  x  10“^ 
coulomb  in  the  MRS  system.  » 

Substitution  of  the  preceding  two  equations  in  the 
equation  immediately  above  involving  and  yields  the 
required  fourth  equation,- 

-rn-  V*  £r  3  e  (D-41) 

Referring  to  the  particle  balance  diagram,  Figure  8, 
in  obtaining  the  time  rate  of  change  of  concentration  equa¬ 
tion  for  each  of  the  three  ions  (by  adding  up  the  contri¬ 
butions  illustrated  for  the  physical  mechanisms  described 
above  )  and  restating  the  preceding  equation  gives  for  the 
four  required  equations: 

|4| -*rn,nx+ 

/)s  <4n. 

f  S  (D-42) 


A 

2-  *r>,Kk,  |£|  +  (D.43) 


-Zjfr  *  0tvx  n,  * 

ut 

(D-44) 

v.  ?  .  *  &-('’< e 
*  ’  ** 

(D-45) 

e. 

In  equations  (D*42) ,  (D-43) ,  (D-44),  and  (D-45), 
the  following  definitions  were  utilized: 

n^  =  electron  concentration,  that  is,  number  of 
electrons  per  unit  volume  at  an  arbitrary 
point  in  the  discharge 
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ti2  =  positive  ion  concentration 

n-j  =  negative  ion  concentration 

q^  =  natural  time  rate  of  ionization  due  to 
natural  causes  such  as  cosmic  rays,  et 
cetera 

a  =  Townsend’s  primary  ionization  coefficient, 
that  is,  the  number  of  ion  pairs  formed  per 
unit  path  length  of  travel  of  an  ionizing 
particle. (electron) 

ar  =  recombination  coefficient  of  electrons  and 
positive  ions 

D^,D2»D3  =  diffusion  constants  for  each  ion 

X  ~  time  for  attachment  of  electrons  to  neutral 
oxygen  molecules  to  form  negative  ions 

=  mobilities  for  each  ion,,  that  is,  ratio  of 
acquired  velocity  to  electrical  field 

E^,E2,E^.  =  composite  field  or  force  per  unit  charge 
acting  on  each  ion 

lSp  =  that  component  of  E  caused  by  electrical 
charge  distribution 

Ew  =  that  component  of  E  caused  by  airstream  or 
*  downwash 

e  =  charge  of  electron 

•Jfc 

dS  =  element  of  area  of  corona  electrode  tip 

i  =  corona  current 

Boundary  conditions  for  the  preceding  four  equations 
are  that  the  solution  for  be  such  that  the  electric  po¬ 
tential  be  constant  over  the  geometry  of  the  corona  electrode 
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(the  surface  of  an  electrical  conductor  an  equipotential) 
and  equal  to  the  applied  voltage.  Also,  Ep  must  be  equal 
to  zero  at  infinity  as  must  np,  n2,  and  n^.  A  similar  but 
less  general  (for  only  two  types  of  ion  and  one  spatial 
dimension)  set  of  equations  is  given  in  Reference  6,  page 
105,  in "the  centimeter-gram- second  system. 


In  theory  (again  computer  methods  would  be  needed 
for  practical  results^),  the  solution  of  the  four  equations 
for  np,  n2,  n^,  and  &p  would  make  possible  determination 
of  the  value  of  corona  discharge  current  resulting  from  an 
application  of  the  designated  magnitude  of  voltage  to  the 
specified  corona  electrode  under  the  given  airstream  con¬ 
ditions.  Jhe  basic  relationship  for  making  this  calcula¬ 
tion  is 


(D-46) 


where  i  is  corona  discharge  current. 


By  definition 


Since  there  are  three  types  of  charge  passing  through 
the  surface  of  the  corona  electrode  (the  electrons,  positive 
ions,  negative  ions),  the  preceding  equation  becomes  a 
vector  sum,  each  term  dealing  with  one  type  of  ion 

*  •  %  *A  %  -  />,  *•  (D'47) 

where  is  volume,  charge  density  for  electrons,  and 

v“p  is  velocity  of  electrons,  et  cetera. 

The  first  and  third  terms  have  minus  signs  because  current 
is  defined  positive  when  in  the  direction  of  positive  charge 
motion.  Using  the  relationships 

y°x  =  nie 

-  n?e 

y°  3  =  n3e 
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and 


vi  =  kiEi 
V2  =  ^2E2 

▼*3  =  k3E^ 

There  follows  that 

effj.-n&X*ni/<xTi-/)yAjTj)  (d-4») 

Evaluation  of  this  expression  by  computer  yields  corona  dis¬ 
charge  current  magnitude. 


As  mentioned  before,  the  above  equations  are  very 
difficult  to  solve.  Hence,  empirical  equations  have  been 
utilized  for  obtaining  the  corona  current.  One  such  em¬ 
pirical  equation  for  corona  discharge  current  given  in 
Reference  2,  page  32  is 

X  =a(v+u/)(v«-vc)  (d-49) 

where  is  applied  voltage  at  corona  electrode, 

Vc  is  corona  threshold  potential, 

A  is  constant, 


v  is  ion  velocity, 

W  is  wind  velocity,  and 

I  is  average  corona  discharge  current. 


The  constant  A  depends  on  ambient  pressure,  tem¬ 
perature,  radius  of  curvature  of  corona  electrode,  and 
point-plane  separation.  No  theoretical  method  has  been 
found  to  determine  the  relationship  of  A  and  the  above 
parameters  and,  hence,  A  must  be  determined  from  experi¬ 
mental  data.  The  velocity  of  ions  v  is  proportional  to 
the  electric  field  strength  and,  consequently,  to  the 
voltage  applied  to  the  electrode.  Using  equation  (D-49), 
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(D-50) 


This  causes  equation  (D-49)  to  become 

I*  a(cv**  w) (v*.-  Vc) 

Thus,  the  corona  current  voltage  characteristics  are  of 
second  degree  for  a  given  wind  speed.  Figure  9  shows 
typical  experimental  values  of  corona  discharge  current 
versus  electrode  voltage  at  two  airstream  velocities  pre¬ 
sented  in  Reference  5  and  from  additional  flight  test  data 
obtained  under  Contract  M  44-177-AMC-114(T) .  The  theoreti¬ 
cal  curve  was  obtained  with  A  =  3.54  x  lO"1^  c  =  2.35  x  10~3, 
and  Vc  =  2000  volts.  It  can  be  seen  that  good  correlation 
exists  between  the  limited  experimental  data  available  and 
theory. 


Experimental  Data  from  Reference  5. 

Probe  Point  Radius  0.002  in 
Corona  Point  Length  24  in 
O  W  =  400  ft/sec 
□  W  =  70  ft/sec 

Unpublished  Experimental  Data  Obtained  under 
Contract  DA  44-177-AMC-114(T) 

AW  =  70  ft/sec 


Figure  9.  Corona  Discharge  Current  Versus  Electrode 
Voltage  for  Two  Wind  Speeds  in  a  Typical 
Discharge  System. 
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SECTION  E:  EFFECTS  OF  ATMOSPHERIC  CONDITIONS 

PRESSURE 


The  operating  pressure  for  the  corona  probe  equals 
the  sum  of  the  air stream  dynamic  pressure  and  the  ambient 
atmospheric  pressure.  The  pressure  level  is  very  important 
to  system  operation  since  the  controlling  field  quantity  at 
constant  temperature  is  E/p  as  described  in  Section  C.  The 
dependence  on  E/p  indicates  that  if  pressure  decreases  by 
some  factor,  the  field,  E,  necessary  to  maintain  the  same 
discharge  current  will  be  decreased  by  the  same  factor  since 
an  accordingly  reduced  E  will  keep  the  value  of  E/p  constant. 
This  is  true  only  up  to  a  certain  limit  of  pressure  reduction 
since  corona  discharge  cannot  occur  in  a  vacuum.  If  the 
value  of  E/p  is  changed  by  variation  of  p,  the  first  Townsend 
coefficient,  described  in  Section  C,  will  go  through  a  maxi¬ 
mum  at  E/p  =  365  volts  per  centimeter  per  millimeter  Hg 
(Reference  6).  Thus,  there  is  a  pressure  for  a  given  corona 
voltage  and  configuration  at  which  corona  current  is  a  maximum. 
Some  control  of  the  airstream  component  of  pressure  to  obtain 
a  possible  optimization  of  performance  is  possible  by  properly 
positioning  the  corona  point  in  the  airstream.  Figure  10 
shows  the  sensitivity  of  discharge  current  to  change  of 
pressure  due  to  altitude. 

TEMPERATURE 


The  effects  of  temperature  can  be  inferred  from  the 
gas  theory  equation. 


P 

N  = 


(E-l) 


where  p  =  pressure, 

N  =  number  of  molecules  per  unit  volume, 
kg  =  Boltzmann’s  gas  constant,  a^d 
T  =  absolute  temperature . 
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Corona  Current,  Microamperes 


II  • 


Figure  10.  Effect  of  Pressure  Change  Due  to  Altitude 
on  Corona  Electrode  Performance. 


This  equation  implies  that  as  temperature  increases  N 
decreases  and  therefore  the  number  of  molecules  available 
for  ionization  decreases.  This  causes  a  loss  in  corona 
current.  If  temperature  is  increased  substantially  (4000°K), 
thermal  ionization  sets  in  and  corona  current  may  begin  to 
increase  with  temperature  (Reference  6). 

MOISTURE 


The  presence  of  humidity  in  the  air  has  a  double 
effect  on  discharger  performance.  It  prevents  the  aircraft 
from  acquiring  charge.  Also,  it  inhibits  the  performance 
of  the  discharger  electrode.  Both  these  effects  are  due  to 
the  tendency  of  water  molecules  to  capture  charge.  Thus, 
in  humid  regions,  the  collection  of  charge  of  both  polari¬ 
ties  into  \rfiich  an  aircraft  flies  discharges  it  over  its 
entire  surface.  Also*  the  operation  of  a  corona  probe  in  a 
high  humidity  atmosphere  results  in  the  formation  of  a 
moisture  layer  over  the  surface  of  the  isolating  supporting 
structure.  This  moisture  layer  presents  an  electrical  path 
between  the  probe  tip  and  the  aircraft  fuselage  which  results 
in  ah  increased  load  of  the  electrical  generator  which  main¬ 
tains  the  probe  potential.  These  generators  are  very  high 
Impedance  devices  and,  hence,  their  terminal  voltage  drops 
due  to  the  increased  load.  This  voltage  drop  results  in  a 
decreased  performance  of  the  discharging  systems. 

PURITY 


Contamination  of  the  corona  point  by  oil,  grease,  or 
tarnish  will  diminish  corona  current  by  reducing  the  second 
ionization  coefficient  for  the  point,  that  is,  secondary 
emission  is  reduced. 
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SECTION  F:  INVESTIGATION  OF  CORONA  PROBE  PERFORMANCE 

IMPROVEMENT 


« 


Reduction  of  corona  electrode  voltage  requirements 
would  be  a  very  significant  contribution  to  electrostatic 
discharger  system  optimization.  Reduced  operating  voltage 
would  permit  a  more  easily  maintained,  smaller,  and  lighter 
discharger  system.  Maximizing  the  magnitude  of  corona 
discharge  current  for  a  given  corona  electrode  voltage 
makes  possible  the  desired  decrease  in  required  operating 
voltage.  This  section  Investigates  several  possible 
approaches  to  the  maximization  of  corona  current. 

THE  RADIOACTIVE  CORONA  PROBE 

Radioactive  materials  emit  particles  which  ionize 
the  surrounding  air,  making  it  conductive.  These  materials, 
therefore,  have  the  ability  to  cause  a  certain  amount  of 
discharge  if  applied  to  an  electrostatically  charged  aircraft. 
However,  the  amount  of  radioactive  material  required  to 
achieve  electrostatic  discharge  if  radioactive  materials  alone 
are  used  is  prohibitively  large  (Reference  2).  The  use  of 
radioactive  material  as  a  method  of  improving  the  performance 
of  the  electrostatic  discharger  system  corona  electrode  is 
conceivable.  The  possible  in^rovements  obtainable  by  use  of 
radioactive  probes  operating  at  high  voltages  can  best  be 
obtained  by  appropriate  experiments. 

IRRADIATION  OF  CORONA  ELECTRODE  WITH  ULTRAVIOLET  LIGHT 

Both  this  method  and  the  use  of  radioactivity  are 
based  on  the  generation  of  an  increased  number  of  triggering 
electrons  for  the  corona  discharge.  Without  some  special 
provision  such  as  these,  the  magnitude  of  corona  current  is 
limited  by  the  relatively  small  ionization  level  set  by 
natural  causes  such  as  cosmic  rays  and  radioactive  traces. 
Ultraviolet  radiation  will  emit  photoelectrons  from  normal 
corona  electrode  materials  such  as  steel  (essentially  iron) 
since  the  energy,  Eu.v.>  of  a  photon  of  ultraviolet  light 
exceeds  the  work- function,  Wf,  of  iron  for  photoelectric 
emission.  The  following  calculation  is  the  basis  for  this 
statement. 
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Eu.v. 

=  6.624  x  10“27  erg-sec  x  10^  cycles/sec 
=  6.624  x  10”^  erg 

=4.16  electron  volts  (F-l) 

where  hp  is  planck  constant,  6.624  x  10-27  erg-sec, 

f  is  ultraviolet  light  frequency,  10^  cps,  and 
Eu#v#  is  energy  of  ultraviolet  light  photon. 

Reference  23,  page  2656,  indicates  that  Wf  for  iron  is 

Wf  -  3.9  electron  volts.  (F-2) 


The  above  calculations:  indicate  that  the  ultraviolet 
photon  energy  is  barely  sufficient  to  produce  electron 
emission  from  iron.  Considering  the  losses  to  be  expected 
due  to  impurities  in  the  probe,  et  cetera,  it  is,  believed 
that  the  use  of  ultraviolet  light  for  increasing  corona  probe 
performance  is  impractical.  •  •• 

COATING  OF  CORONA  ELECTRODE  WITH  PHOrOEMISSIVE  MATERIAL 

If  the  corona  electrodes  -are  coated  wit*1  °>atarlals 
especially  selected  for  low  photoelectric  work  function  such 
a 2  activated  cesium,  Irradiation  of  the  corona  electrode  with 
light  in  the  visible  range  could  result  in  -performance  improve¬ 
ment.  Suitable  experiments  should  be  performed  to  determine 
the  extent  of  these  improvements. 

OPERATION  OF  CORONA  ELECTRODE  IN  ENGINE  EXHAUST 

In  principle,  an  engine  exhaust  is  a  possible  means 
of  accomplishing  alone  or  assisting  a  corona  discharge  system 
in  removal  of  electrostatic  charge  from  aircraft.  This  is 
due  to  the  presence  of  ions  in  the  hot  gasses  of  the  engine 
exhaust  caused  by  the  phenomenon  of  thermal  Ionization. 

These  ions  improve  the  performance  of  the  corona  electrode  by 
contribution  to  the  number  of  trigger  electrons  in  a  way  sim¬ 
ilar  to  the  mechanism  of  current  enhancement  by  radioactivity. 
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Unfortunately,  the  practical  use  of  engine  exhaust 
Ionization  to  ln^prove  electrostatic  discharger  system  per¬ 
formance  Is  very  severely  limited  by  the  exhaust  temperature 
requirements.  The  results  of  calculation  using  the  basic 
equation  of  thermal  Ionization,  Saha's  equation  (Reference  6), 
verify  the  excessive  temperature  requirements.  Saha's 
equation  for  one  component  of  the  engine  exhaust  Is 

f  *  *te*,o-Tr™  e  *  <F-3) 

/-rx 

where  p  Is  total  pressure  In  atmospheres, 

T  Is  the  temperature  of  the  exhaust  In  degrees, 
kelvln, 

eVl  Is  the  energy  required  for  Ionization  of  the 
gas  atoms  In  ergs, 

kg  Is  Botzman's  constant,  and 

r  Is  fraction  of  atoms  Ionized  In  each  mole  of 
engine  exhaust  as  It  leaves  engine. 

We  calculate  the  value  of  r  for  CO2,  a  typical 
engine  exhaust  component  (Reference  6)  at  the  following  con¬ 
ditions  : 

P  -  1  atmosphere 

T  =  2000°K  (3140°F) 

Vt  =  5.5  volts. 

Substituting  values  In  Saha's  equation, 


therefore. 


r  0 
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The  low  value  of  r  obtained  even-  at  elevated  tem¬ 
perature  demonstrates  the  nonfeasibility  of  corona  probe 
performance  improvement  by  the  use  of  helicopter  engine 
exhaust.  The  operating  temperature  of  helicopter  engine 
exhausts  (Reference  21)  is  approximately  800°F.  The  tempera¬ 
ture  of  jet  engine  exhaust  is  also  too  low  for  use  in  electro¬ 
static  discharger  system  performance  improvement  with  jet 
aircraft . 

CHOOSING  CORONA  ELECTRODE  MATERIAL  FOR  SECONDARY  EMISSION 
CAPABILITY 


The  role  of  secondary  emission  phenomena  in  the 
generation  of  corona  current  is  represented  by  the  parameter, 
y ,  discussed  in  Section  C,  page  13*  As  stated  there,  the 
effect  of  secondary  emission  in  the  physical  process  of  corona 
discharge  does  not  become  significant  until  the  latter  stages 
of  the  air  electrical  breakdown*  Since  electrical  breakdown 
in  a  corona  discharge  for  a  discharging  system  must  be  main¬ 
tained  far  below  the  electrical  breakdown  condition  (spark), 
the  use  of  secondary  emission  techniques  to  enhance  corona 
probe  performance  is  not  considered  practical. 

CORONA  POINT  GEOMETRY 

As  has  been  previously  discussed,  corona  point 
geometry  and  location  on  the  fuselage  significantly  affect 
corona  point  performance.  The  location  of  the  probe  in  the 
maximum  downwash  velocity  region,  as  suggested  in  Figure  1, 
markedly  Increases  corona  discharging  current.  Minimum 
corona  point  radius  of  curvature,  as  discussed  in  Section  C, 
also  Increases  corona  performance.  Practical  design  aspects, 
however,  limit  this  radius  of  curvature  within  0.002  to 
0.006  inch. 

The  effect  of  corona  point  geometries  on  discharging 
current  performance  has  also  been  investigated  in  Reference 
22.  Figure  11,  reproduced  from  Reference  22,  compares  the 
performance  of  several  different  corona  point  configurations. 
As  can  be  seen  from  this  figure,  sharp  needle  probes  Increase 
the  probe  performance.  The  use  of  such  points,  however,  is 
impractical  in  operational  aircraft  because  of  possible 
injury  to  ground  personnel. 
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Discharge  Current,  I^*a 


Probe  Type 

Wick 

Sharp  Needle  Probe 
Hollow  Cupped  Needle  Probe 
Solid  Cupped  Needle  Probe 
Radial  Cluster  Sharp  Needle  Probe 
Radial  Cluster  Hollow  Needle  Probe 

Knife  Edge  Probe 

1/4-in-Diameter  Cupped  Probe 
Thin  Wire  Probe 

Figure  11.  Corona  Current  Discharge  for  Various 
Corona  Point  Configurations  (60-Kllo- 
volt  Applied  Voltage  and  90-railes-per 
hour  Velocity). 
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SECTION  G:  DISCHARGER  SYSTEM  OPERATION  AND  ANALYSIS 


THE  ELECTROSTATIC  DISCHARGER  PRINCIPLE 

In  the  absence  of  a  discharger  system,  triboelectric 
or  other  charging  mechanisms  will  cause  the  electrostatic 
charge  and,  consequently,  the  potential  of  the  aircraft  with 
respect  to  its  environment  to  rise  steadily.  When  the  result¬ 
ing  electrical  field  strength  in  the  atmosphere  adjacent  to 
the  aircraft  reaches  a  critical  value  of  approximately  30 
kilovolts  per  centimeter,  natural  corona  current  from  random 
points  of  the  aircraft  surface  occurs.  This  corona  current 
is  the  flow  of  charge  generated  by  ionization  of  the  atmos¬ 
phere  which  occurs  when  the  electrical  field  becomes  suffi¬ 
ciently  large,  that  is,  reaches  the  critical  value  (see  Section 
C).  The  rise  of  aircraft  potential  is  halted  only  when  this 
flow  of  charge  off  the  airframe  equals  the  Incoming  natural 
charging  current.  The  value  of  aircraft  potential  at  which 
corona  occurs  is  determined  by  the  sharpness,  that  is,  radius 
of  curvature  of  the  points  on  the  aircraft  surface  and  there¬ 
fore  varies  with  the  aircraft.  A  5-kilcv^lt  potential  is  a 
typical  threshold  voltage,  causing  corona  discharge  from 
points  of  0.06-inch  radius  of  curvature.  The  stray  discharge 
current  increases  as  the  triboelectric  aircraft  voltage  moves 
further  away  from  zero  until  the  total  corona  discharge 
current  equals  the  charging  current,  that  is,  until  the  net 
charging  current  is  zero.  The  amount  of  charge  on  the 
aircraft  becomes  constant  at  this  time.  There  is  a  simple 
relationship  between  the  charge  on  the  aircraft  and  its 
potential  provided  that  the  earth  is  considered  the  only 
other  conductor  present  in  addition  to  the  aircraft.  When 
this  quite  valid  approximation  is  made 

VH  =  Qh  (G-l) 

where  VH  =  aircraft  potential  relative  to  earth, 

Qh  *  aircraft  charge,  and 
Cjj  *  aircraft  capacitance. 
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The  value  of  Ch  varies  inversely  with  the  altitude  and 
directly  with  the  area  of  the  aircraft  and  is,  therefore, 
not  constant  during  flight.  Techniques  for  measuring  Ch 
are  given  in  Reference  5.  The  problem  arising  from  air- 
craft  electrostatic  charging  is  due  to  the  magnitude  of 
the  residual  charge  remaining  on  the  aircraft  after  dis¬ 
charge  by  natural  process. 

The  dynamic  discharger  system  described  in  this 
report  prevents  excessive  rise  of  potential  by  auto¬ 
matically  controlling  the  magnitude  of  corona  discharge 
current  from  one  deliberate  corona  discharge  point  on  the 
aircraft.  Figure  12  illustrates  the  system  and  the  charge- 
discharge  current  loop.  The  design  and  construction  of 
this  discharger  electrode  gives  it  a  superior  corona  dis¬ 
charge  capability.  Consequently,  the  system  rapidly  brings 
the  net  charge  current  to  zero  with  small  residual  aircraft 
charge  and  potential.  Corona  occurs  only  at  the  discharger 
probe  because  the  discharger  system  maintains  the  aircraft 
at  a  potential  Substantially  lower  than  the  corona  thresh¬ 
old.  The  control  of  discharge  current  is  accomplished  by 
setting  the  discharge  electrode  at  a  potential  difference 
with  respect  to  the  aircraft  surface  by  means  of  a  variable 
voltage  generator.  This  voltage  difference  determines  the 
field  strength  about  the  electrode  and,  hence,  the  discharge 
current  as  explained  in  Section  D.  In  effect,  the  discharger 
system  ionizes  the  atmosphere  at  the  corona  electrode  creating 
a  conductive  leak- off  path  at  this  point  for  the  charge 
accumulated  on  the  aircraft.  The  ions  of  the  corona  dis¬ 
charge  current  satisfy  the  closed-path  law  of  electric  current 
by  eventually  returning  to  the  particles  (dust,  snow,  et 
cetera)  from  which  they  were  triboelectrlcally  removed  (see 
Figure  12).  Since  natural  charging  currents  of  both  polar¬ 
ities  are  possible,  the  system  must  be  capable  of  supplying 
both  positive  and  negative  voltage  to  the  corona  discharge 
electrode.  The  polarity  of  voltage  applied  to  the  corona 
electrode  Is  that  of  the  electrostatic  charge  acquired  by 
the  aircraft.  Thus,  both  positive  and  negative  corona 
characteristics  are  pertinent  to  the  discharger  operation. 
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Triboelectrically 
Charged  Particles 


•  • 


12.  Block  Diagram,  Aircraft  Static  Electricity  Discharging  System. 


DISCHARGER  SYSTEM  OPERATION 

A  fully  effective  discharger  system  has  been  con¬ 
structed  by  Dynasciences  Corporation  under  USAAML  Contract 
DA  44-177-AMC-114(T) .  It  is  basically  a  feedback  amplifier 
using  a  sensor  of  the  generating  voltmeter  type  (see  system 
diagram,  Figure  13).  In  operation,  aircraft  charge  is 
acquired  at  the  rate  In  predominantly  by  triboelectric 
effects  on  the  moving  rotor  blades  of  th?  aircraft.  This 
causes  the  aircraft  to  assume  a  potent ir vr,  with  respect 
to  its  environment.  In  most  cases,  the  earth’s  electrostatic 
field  is  small  enough  to  consider  this  environmental  poten¬ 
tial  as  equal  to  ground  potential.  The  acquired  charge  also 
causes  an  electrical  field,  Er,  which  is  a  measure  of  the 
potential,  ^H,  to  exist  about  the  aircraft.  The  generating 
voltmeter  sensor  develops  an  A.C.  signal,  vs>  proportional 
to  EH  (and  hence  to  VH  as  shown  on  page  85  )  which  is  the 
input  to  the  control  unit.  This  unit  essentially  establishes 
the  magnitude  and  polarity  of  the  corona  electrode  voltage 
by  the  magnitude  of  its  output,  Vc,  which  will  increase  or 
decrease  with  In.  The  control  unit  also  contains  a  compen¬ 
sator  circuit  which  Imposes  the  correct  transient  response. 
Without  the  compensator  circuit,  the  transient  response 
would  be  determined  by  the  relatively  large  time  constant, 
tg,  of  the  high  voltage  generator.  The  system  would  oscillate 
excessively  or  "hunt"  if  not  for  this  corrective  circuit. 

The  power  amplifier  and  high  voltage  generator  develop  the 
large  D.C.  corona  electrode  voltage  (in  the  kilovolt  range) 
corresponding  to  vc  and  the  charging  current,  In.  This  large 
D.C.  potential  causes  the  atmosphere  about  the  discharger 
probe  to  ionize  and,  consequently,  the  acquired  charge  leaves 
the  aircraft  at  this  point.  The  effective  charging  current, 
therefore,  is  In  -  Id  =  It*  The  system  automatically  increases 
the  corona  electrode  voltage  until  Id  =  In  or  It  -  0.  The 
aircraft  electrostatic  charge  ceases  to  rise  at  the  time  the 
net  charging  current,  Itf  becomes  zero  since  charge  Is  then 
leaving  the  aircraft  as  fast  as  it  is  impinging  on  it.  The 
amount  of  charge  on  the  aircraft  remains  constant  once  It 
becomes  zero.  In  combination  with  the  aircraft  capacitance, 

Cr,  this  residual  charge  determines  the  magnitude  of  aircraft 
potential,  Vr.  Proper  system  design  using  the  results  of  the 
loop  analysis  on  page  78  insures  that  the  final  value  of  air¬ 
craft  potential  corresponds  to  an  energy  level  below  the 
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danger  point  (one  mill! joule)  for  the  given  charging  con¬ 
ditions.  This  condition  stated  mathematically  is 

£a  *  /x  Jo  vie  (g- 

where  EA  =  electrostatic  stored  energy  of  aircraft, 

Chm  =  maxiraum  aircraft  capacitance  (at  lowest 
operating  altitude),  and 

VHF  -  final  aircraft  potential. 

The  value  of  is  set  by  the.  system  gain  as  shown  in  the 
loop  analysis.  The  discharger  is  a  servo  system;  because 
its  output,  I^,  affects  the  magnitude  of  its  effective 
input,  It. 


The  magnitude  of  corona  discharge  Increases  with 
the  velocity  of  the  air stream  sweeping  past  the  corona 
electrode  as  discussed  in  Section  D.  Consequently,  this 
electrode  is  positioned  under  the  rotors  so  that  the  wind 
of  the  rotor  downwash  improves  probe  performance  as  much 
as  possible.  Figure  14  on  page  76  shows  the  mounting  of 
the  electrode  on  the  fuselage  of  a  tandem  rotor  helicopter. 

The  operating  effectiveness  of  the  electrostatic 
discharger  system  can  be  checked  while  in  flight  at  the 
pilots  cockpit  and  at  the  cargo  area  by  means  of  the 
failsafe  feature.  This  consists  of  two  control  units 
each  having  red  and  green  indicator  lights.  One  unit  is 
in  the  cockpit,  the  other  is  at  the  cargo  area.  When 
the  discharger  system  is  satisfactorily  preventing  the 
rise  of  aircraft  electrostatic  charge,  the  green  light  is 
on  in  both  control  units.  The  red  light  will  go  on  if  the 
sensor  is  defective,  that  is,  if  the  sensor  voltage  output 
is  zero.  The  red  light  will  also  go  on  if  the  sensor  out¬ 
put  is  above  a  specified  voltage  level.  Such  an  over¬ 
voltage  indicates  an  excessive  amount  of  electrostatic 
charge  on  the  aircraft  because  of  failure  of  some  discharger 
system  unit  or  extreme  operating  conditions.  When  the 
system  is  Initially  energized  in  charging  conditions,  the 
red  light  will  go  on.  However,  when  the  system  has  dissi¬ 
pated  the  excess  charge  built  up  before  the  system  was 
switched  on,  the  green  light  will  normally  come  on. 


•CORONA  DISCHARGE 
ELECTRODE 


Figure  14.  Typical  Mounting  of  Corona  Discharge 

Electrode  on  a  Tandem  Rotor  Helicopter 


LOOP  ANALYSIS  OF  DISCHARGER  SYSTEM 


The  system  analyzed  here  is  the  generating  voltmeter 
sensor  system  (GVSS).  The  purpose  of  the  analysis  is  to  find 
the  time  domain  system  response,  Vft,  to  a  step  input,  In.  To 
do  this,  we  first  find  the  system  parameter,  V){/In,  in  the  s 
domain  and  then  apply  the  Inverse  laplace  transform  to  this 
parameter  to  find  the  time  response.  The  required  system 
parameter  is  written  as  the  product  of  two  other  more  easily 
calculated  system  parameters.  The  analysis  proceeds  as 
follows,  with  reference  to  Figure  13. 

—  -  Xc  (G-3) 

The  first  parameter  evaluated  is  I^/In*  Accordingly, 

Zt  *  -T*  -  Zp  (G-4) 


X*  I*  / 

Xp  x++Zp  +  / 


(G-5) 


The  system  parameter,  appearing  in  equation  (G-5)  is 

the  open  loop  gain,  Fol>  given  by 


K _ 

C0%  ft  *  -6cs  ) 


(G-6) 


where 


(G-7) 


Applying  the  expression  obtained  for  Id/It  t0  equation  (G-4), 
one  obtains  the  following  result  for  Id/In,  the  closed  loop 
gain; 

r  / _ 

+  /  (G-8) 

X 
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(G-9) 


<v' 


X _ 

CMS(/+tcs)  +  K 


This  evaluates  one  of  the  two  required  system  parameters. 
The  remaining  one  Is  obtained  by  taking  the  reciprocal  of 
Id/VH,  as  follows. 


Z/>  _  X 

*  t  +  *cS 


(G-10) 


Thus,  substituting  in  equation  (G-2) 

_  / +  *CS  ,  X _ 

JC/f  X  (/  +  ic*)  ♦  n 


(  X  *  X 


(G-ll) 

(G-12) 


This  is  the  required  network  parameter  in  the  s  domain. 

The  laplace  transform  of  the  system  response  to  a 
step  input  of  In  is  obtained  by  multiplying  VH/In  by  the 
current  input  transform  K^/s,  where  is  the  step  amplitude. 
Accordingly,  the  system  response,  Vjj,  to  an  input  step  of 
In  is 

vy.  **0**c*) 

l  C*^S*  f  - *  /  J  (0-13) 


Using  Reference  13,  the  time  domain  expression  for  Vh/Kx 
derived  from  this  is 


Kx 


e. 


14) 
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where 


■0"  *  |u/  J/m  {* J*  + 


f  *  w  1*  f  wlnAL 


This  is  the  required  result,  that  is  the  system 
response,  Vh,  to  a  step  input  of  charging  current,  In. 

From  the  solution,  it  is  seen  that  in  the  steady  state, 
that  is  when  t  becomes  infinite 

ft* 

*  <:  (G-15) 

This  relationship  sets  the  value  of  overall  system  gain, 

K,  for  that  electrostatic  discharger  system  idiich  brings 
the  aircraft  potential  to  a  specified  value  for  a  given 
charging  current. 

Also  of  interest  is  the  value  of  discharge  current, 
Id,  in  the  steady  state  after  a  step  input  of  natural 
charging  current,  In,  has  been  applied  to  the  system.  This 
is  obtained  by  the  use  of  the  expression  for  closed  loop 
gain,  F„  given  by  equation  (G-8).  Rewriting  this 
expression 

_  b  9  k 

c ’>•  X „  C^S(T**7sTTK  (G-16) 


The  laplace  transform  of  1^  when  a  step  input  of  In  is 
applied  to  the  system  is  obtained  by  multiplying  Id/In  by 
the  appropriate  current  input  transform,  K^/s.  K*  is  the 
step  amplitude.  Accordingly,  the  laplace  transform  of  1^ 
to  such  an  input  is 
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X,.  -**» _ 

*  k] 


(G-17) 


The  time  domain  expression  for  I<j  derived  from  equation 
(G-16)  is  given  by  Reference  13  as 


X0  •  KA 


/  -  _  "I 

'* pi* e  ■*<e  «*(# /'-  f  *-  fyl 


(G-18) 


where 


(J  * 


J£*c 

V 


•i 


The  fact  that  I&  becomes  equal  to  the  magnitude  of  charging 
current,  K~,  in  the  steady  state,  that  is  when  t  ■  •€> 
verifies  tne  effectiveness  of  the  system  as  an  electrostatic 
discharger. 

An  analog  computer  analysis  of  the  discharger  system 
dynamics  is  contained  in  Appendix  II. 


DESIGN  CURVES  FOR  ELECTROSTATIC  DISCHARGER  SYSTEM 

The  quantities  necessary  to  completely  specify  a 
discharger  system  design  are: 

1.  System  steady  state  gain,  K,  that  is, 
ratio  of  charging  current  to  steady 
state  aircraft  potential  Introduced 
on  page  77. 
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2.  Compensator  components. 

3.  Maximum  corona  electrode  voltage  which 
system  can  apply. 

4.  Number  of  corona  electrodes  required. 

The  magnitude  of  these  quantities  in  any  particular  dis¬ 
charger  design  is  set  by  the  existing  values  for  the 
following  operating  conditions: 

1.  Maximum  charging  current. 

2.  Maximum  aircraft  potential  without  corona. 

3i  Aircraft  capacitance  at  operating  altitude. 

4.  Aircraft  speed. 

5.  Time  constant  of  high  voltage  supply^. 

6.  Required  damping  ratio. 

Curves  enabling  quick  determination  of  the  system 
specifications  for  given  charging  conditions  are  given  in 
Figure  15.  The  curves  give  the  required  steady  state  gain 
K,  compensator  components,  Rl,  R2,  and  Cc,  and  corona  elec¬ 
trode  voltage,  Vp,  over  a  ranee  of  charging  conditions  for 
specified  damping  ratio  (0.35)  and  aircraft  electrostatic 
energy  (one  mllll joule).  The  curves  are  entered  with  the 
following  quantities:  natural  charging  current;  In,  wind 
speed;  W,  helicopter  capacitance;  Ch,  and  high  voltage 
generator  time  constant;  t».  Appendix  I  shows  how  the 
data  for  the  curves  are  calculated.  The  units  to  be  used 
with  the  curves  are: 

In  -  microamps 

C|{  -  micromicrofarads 

K  -  microaraps  per  kilovolt 

W  -  feet  per  second 

Vp  -  kilovolts 
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seconds 


tc  -  seconds 
R2  -  arbitrary 

The  curves  are  used  as  follows. 


1.  Enter  with  In. 

2.  Read  Vp  using  appropriate  wind  speed  curve. 

3.  Read  K  using  appropriate  solid  helicopter 
capacitance  curve  In  upper  left  comer. 

4.  To  determine  Rl/R2?  use  appropriate  dashed 
helicopter  capacitance  curve  (evaluates 

In/  Ch)  in  upper  left  corner.  Then  proceed 
to  the  appropriate  t»  curve  in  the  upper 
right  corner.  Read  R1/R2  at  the  inter¬ 
section  of  the  vertical  at  this  point  with 
the  R1/R2  axis. 

5.  To  determine  CcR2,  project  the  vertical 
established  in  step  4  downward  to  Intersect 
the  appropriate  tg  curve  in  the  lower  right 
hand  corner.  Draw  a  horizontal  at  this 
point  of  intersection.  Read  CcR2  at  the 
intersection  of  the  horizontal  with  the 
CcR2  axis. 

6.  Multiply  R1/R2  and  divide  CcR2  by  the 
selected  value  of  R2  to  obtain  R1  and  Cc. 


SENSORS  IN  ELECTROSTATIC  DISCHARGER  SYSTEMS 


The  function  of  the  sensor  is  to  provide  electrical 
Information  to  the  discharger  system  concerning  the  magni¬ 
tude  of  the  aircraft  potential.  The  theory  and  relative 
merits  of  three  sensors  are  discussed  below.  The  first  of 
these,  the  generating  voltmeter  sensor,  has  been  used  in 
the  discharger  system  because  of  its  superior  qualities. 
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A  rotating  vane  generating  voltmeter  produces  an 
A.C.  output  voltage  with  an  amplitude  directly  proportional 
to  the  field  strength  of  the  charged  body  to  which  it  is 
attached.  When  attached  to  the  surface  of  an  electrostat¬ 
ically  charged  aircraft,  this  device  senses  the  aircraft 
potential  since,  as  proven  below,  the  field  strength  is 
directly  proportional  to  the  aircraft  potential.  This 
sensor,  however,  is  ineffective  at  a  point  on  the  aircraft 
which  is  in  corona.  At  such  a  point,  charge  density  no 
longer  builds  up  as  aircraft  potential  rises;  instead,  it 
leaks  off  into  the  atmosphere.  Hence,  the  electric  field 
at  such  points  is  not  proportional  to  aircraft  potential. 
This  type  of  sensor  should  therefore  be  located  on  the 
aircraft  at  a  point  of  minimum  curvature.  The  surface  must 
be  convex  at  the  point  of  location,  that  is,  the  curvature 
must  be  positive  in  order  for  the  field  strength  to  be 
strong  enough  to  cause  a  sensor  response.  When  located  at 
such  a  point,  the  output  of  this  sensor  is 

Kr  s  Kra*  (G-19) 

where  *  field  of  aircraft  surface  at  sensor 

location  and 

K$en  *  transfer  function  of  sensor. 

It  is  now  shown  that  if  the  sensor  is  placed  at  a 
point  on  the  aircraft  which  is  roughly  spherical  with  a 
radius  of  curvature  of  1  meter,  the  constant  of  proportion¬ 
ality,  Ksen>  relating  the  field  strength  to  the  aircraft 
potential  has  a  magnitude  very  near  unity.  To  see  this, 
consider  the  potential  at  such  a  location  on  the  aircraft 
surface.  Assuming  that  the  surface  can  here  be  approxi¬ 
mated  by  an  electrically  Isolated  charged  sphere, 

V  .  J_ _ _ 

*  r/re0At  (G-20) 

where  q  =  charge  on  sphere, 

£0  =  permittivity  of  free  space,  (MKS  units). 
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Rs  =  radius  of  curvature,  and 


VH  =  aircraft  potential. 


Multiplying  this  ratio  by  Rs/Rs  gives 


VW 


q  *x 


(G-21) 


According  to  Reference  8,  the  electric  field  at  the 
surface  of  an  aircraft  approximated  by  an  isolated  sphere 
of  radius,  Rg,  is 


_ 

if/re9  a/ 


Therefore,  equation  (G-20)  for  Vg  becomes 


(G-22) 


v*  ‘  At 


(G-23) 


where 


Thus,  if  Rs 


=  field  strength  at  that  point  on  the 
aircraft  where  the  curvature  is 
approximately  constant  with  radius,  Rs. 

=  1  meter. 


Yt  •  A  ft#  8  Affff  .y.ii  m  Kfgff  Vg  (G- 24) 

*x 

Passive  Corona  Point.  Sensor 

The  system  employing  this  type  of  sensor  is  the  same 
as  the  generating  voltmeter  sensor  system  described  on  page  84 
except  that  the  amount  of  corona  current  passing  through 
a  passive  corona  point,  that  is,  one  connected  directly  to 
the  aircraft,  serves  as  a  sensor  of  the  aircraft  potential. 

The  system  diagram  is  that  of  the  GVSS  system  shown  on  page 
87  except  for  the  sensor  transfer  function. 
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In  order  for  the  passive  corona  point  to  go  into 
corona,  the  aircraft  potential  roust  typically  rise  to  a 
potential  of  3  to  5  kilovolts.  Figure  16  is  a  plot  of  the 
transfer  function  of  this  type  of  sensor.  Thus,  this  type 
of  system  is  incapable  of  reducing  the  aircraft  potential 
below  3  to  5  kilovolts.  Since  this  corresponds  to  an 
electrostatic  energy  level  above  the  maximum  safe  value  of 
one  millijoule  for  most  operating  conditions  (helicopters 
over  7000-pounds  gross  weight  hovering  below  twenty- five 
feet),  this  system  is  unacceptable  (Reference  4). 


Dynamic  Nentrali zer_Sen&gr 

This  device  is  actually  a  combination  of  a  potential 
sensing  method  and  a  system  concept.  A  loop  analysis  appears 
in  Figure  17. 

In  this  system,  two  high  voltage  generators  are 
connected  with  opposite  polarity  to  the  aircraft  frame  at 
one  terminal  and  to  a  corona  electrode  at  the  other  terminal 
respectively,  as  shown  in  Figure  12  on  page  72  •  The  voltage 
output  of  the  generators,  VcpOS  and  Vcneg,  is  adjusted  until 

~  V€Mt*  /<0gr  *  0  (G- 25) 

where  Kpos  =  transfer  function  of  positive  corona 

point,  and 

Kneg  =  transfer  function  of  negative  corona 
point . 

As  described  in  Section  C,  KpOS  does  not  equal  Kneg 
because  of  differences  in  the  mechanisms  of  positive  and 
negative  corona.  To  prevent  other  differences  in  the  corona 
electrode  operating  conditions,  the  relative  positions  of 
both  corona  electrodes  with  respect  to  the  aircraft  and  the 
air  speed  are  made  identical.  Under  these  circumstances, 
and  with  the  aircraft  potential  zero,  the  discharge  current 
from  both  corona  electrodes  will  be  equal  and  the  net  dis¬ 
charge  current  will  be  zero.  If  the  aircraft  potential 
becomes  positive  (for  the  acquisition  of  negative  potential 
all  polarities  should  be  reversed  in  the  following  dis¬ 
cussion)  the  total  potential  difference  between  the  positive 
corona  electrode  and  the  surrounding  environment  increases. 
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Also,  the  total  potential  difference  of  the  negative  corona 
electrode  with  respect  to  its  environment  decreases.  Con¬ 
sequently,  the  output  of  the  positive  corona  probe  increases 
and  the  output  of  the  negative  corona  probe  decreases.  As  a 
result  of  the  imbalance  of  the  positive  and  negative  currents, 
the  aircraft  takes  in  a  net  charge  with  a  polarity  that  is 
opposed  to  that  of  the  natural  charge  which  created  the  air¬ 
craft  potential  (Reference  4).  The  system  therefore  dis¬ 
charges  the  aircraft. 

A  loop  analysis  of  the  dynamic  neutralizer  system 
follows.  The  mathematical  techniques  used  in  the  analysis 
of  the  GVSS  system  on  page  77  are  also  used  here.  The 
system  diagram  used  is  shown  in  Figure  17. 

The  purpose  of  the  analysis  is  to  find  the  time 
domain  response,  V^,  to  a  step  input  of  In  for  this  system. 
The  required  system  parameter,  VH/In,  is  written  as  the 
product  of  two  other  system  parameters.  Accordingly, 


(G- 26) 


The  first  parameter  evaluated  is  Id/In»  We  use  the  following 
relationship  for  net  charging  current. 

If  -  I/f  “  I D  (G- 27) 

Manipulation  of  equation  (G-27)  yields 


Ip 

In 

To  evaluate  Id/In» 
follows . 


(G-28) 


Id/It  is  obtained  from  the  analysis  which 


Ip  *  I*ot  +  I*te  (G-29) 

s  (y» f  (Vjt  m  Vgm*)  (G-30) 

r  *  X*,c)  +  W  -  W  (G”3I) 
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Xt 

G„S 


r  *  -£*- 
*  c** 


'*>r  T  r'  nr* 


and,  therefore, 


X*  *  tt 

Xt  cMs 


Hence,  substituting  in  equation  (G-28) 


f  A***  *  ***» 


(G-33) 

(G-34) 

(G-35) 

(G-36) 


Prom  equation  (G-32) 


*. 

/ 

(G-37) 

since 

Jk. 

Kr  X* 

— •  »  ■ 

(G-38) 

X,  X„ 

Applying  equations  (G-37)  and  (G-36)  to  equation  (G-38) 
gives 


(G-39) 


The  system  response  to  a  step  input  of  I  is  obtained 
by  multiplying  VH/l  by  the  step  input  transform?  K/s,  where 
is  the  step  amplitude*  Accordingly,  the  system  response, 
Vjj,  to  an  input  step,  In,  is 
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(G-40) 


CjSx+  *  **g€ )  S 


The  time  domain  expression  for  V^/K  derived  from  this  is 


Km  K/mtff  + 

Equation  (G-41)  indicates  that  the  dynamic  neutra¬ 
lizer  is  an  inherently  stable  system.  However,  as  has  been 
shown  in  Reference  4,  the  gain  of  this  system  is  insufficient 
to  maintain  the  helicopter  energy  within  one  milli joule  at 
charging  rates  above  0.5 


/-  e 


5 


V 


(G-41) 
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CONCLUSIONS 


The  following  conclusions  were  reached: 

1.  The  theory  developed  and  included  in  this 
report  predicts  discharge  system  operation 
consistent  with  its  actual  operating  charac¬ 
teristics.  It  therefore  puts  the  discharger 
system  design  on  a  scientific  rather  than 
empirical  basis. 

2.  Analytical  procedures  providing  guidelines 
for  discharger  system  design  and  optimiza¬ 
tion  have  been  established  by  this  program. 

3.  To  prevent  the  loss  of  discharger  system 
effectiveness  caused  by  return  of  discharge 
current  to  the  aircraft,  the  system  dis¬ 
charge  electrode  should  be  mounted  so  that 
the  airstream  past  the  electrode  does  not 
strike  any  part  of  the  aircraft. 

4.  The  amount  of  polarization  in  helicopter 
engine  exhaust  is  too  low  to  improve  dis¬ 
charger  system  performance. 


RECOMMENDATIONS 

Several  approaches  to  discharger  system  performance 
Improvements  discussed  in  Section  F  of  this  report  should  be 
evaluated  experimentally.  Specifically,  it  is  recommended 
that  appropriate  tests  be  performed  to  determine  the  in¬ 
crease  of  corona  point  performance  obtainable  by  the  use  of 
radioactive  or  photoemissive  materials. 
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US  Army  Combat  Developments  Command  Transportation  Agency  1 

US  Army  Combat  Developments  Command  Experimentation  Command  1 
US  Army  War  College  1 

US  Army  Command  and  General  Staff  College  1 

US  Army  Transportation  School  1 

US  Army  Aviation  School  1 

US  Army  Quartermaster  School  1 

US  Army  Transportation  Center  and  Fort  Eustis  1 

US  Army  Infantry  Center  1 

US  Army  Aviation  Test  Board  1 

US  Army  Arctic  Test  Center  1 

US  Army  Electronics  Command  ,1 

US  Army  Aviation  Test  Activity  1 

US  Army  Transportation  Engineering  Agency  1 

Air  Force  Flight  Test  Center,  Edwards  AFB  1 

Air  Proving  Ground  Center,  Eglin  AFB  1 

Air  Force  Avionics  Laboratory,  Wright-Patterson  AFB  1 

Chief  of  Naval  Research,  D/N  1 

US  Naval  Air  Station,  Patuxent  River  1 

Marine  Corps  Educational  Center  1 

Marine  Corps  Liaison  Officer,  US  Army  Transportation  School  1 

Testing  and  Development  Division,  U[S  Coast  Guard  1 

Ames  Research  Center,  NASA  1 

Lewis  Research  Center,  NASA  1 

Manned  Spacecraft  Center,  NASA  1 

NASA  Representative,  Scientific  and  Technical  Information 

Facility  2 

Research  Analysis  Corporation  1 

NAFEC  Library  (FAA)  1 

Defense  Documentation  Center  20 

US  Government  Printing  Office  1 
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The  values  of  K,  Rl,  Cc,  and  Vp  required  to  obtain 
the  curves  are  found  by  simultaneously  solving  the  six 
equations  which  apply  in  terms  of  In,  W,  Cjj,  and  tg.  These 
equations  are 
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Equation  (APPI-4)  is  a  linearization  of  the  discharge 
current  -  wind  speed  dependence  at  a  given  Vp.  Equations 
(APPI-5)  and  (APPI-6)  come  from  the  condition  that  the  com¬ 
pensator  transfer  function  be  given  by 
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for  the  compensator  circuit  as  shown  on  the  curves.  The 
constants  in  the  above  equations  and  their  values  used  in 
constructing  the  curves  are  as  follows. 


€ 
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=  permissible  helicopter 
electrostatic  energy 

=  system  damping  ratio 


10~3  joule 


0.35 


K(jp  transfer  function  of 

corona  point  (a  two-point  0.835  Mamp/Kv 

probe)  ^ 


Wc,  value  of  wind  at  which  80  ft/ sec 

K(jp  was  measured 


Solving  the  equations  and  applying  suitable  con¬ 
stants  to  allow  the  use  of  convenient  units  yields 
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APPENDIX  II 


ANALOG  COMPUTER  ANALYSIS  OF  DISCHARGER  SYSTEM  DYNAMICS 

An  analog  computer  analysis  of  the  discharger  system 
dynamics  was  performed.  The  system  equations  are  as  follows. 
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The  symbols  are  defined  in  the  block  diagram  shown  in  Figure 
18  of  this  appendix.  Values  for  the  system  constants  are  as 
follows . 
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Cjj  =  10  ^  farads 

Kcpu  -  3  X  10"3  VOLTS  PEAK 

VOLTS  D.C. 


(APPII-8) 

(APPII-9) 


K  =  to  be  varied  between  2.67  (APPII-10) 

1  X  10"  3  and  13.31  X  10‘3 
VOLTS  PEAK 
VOLTS  PEAK 

K  =  103  VOLTS  PEAK  (APPIJ-11) 

VOLTS  PEAK 

K  =2.5  X  103  VOLTS  DC  (APPII-12) 

8  VOLTS  PEAK 


tg  =  to  be  varied  between  0.03  (APPII-13) 
and  0.5  second 

Kcp  =  0.835  X  10“9  AMPERES  DC  (APPII-14) 

VOLTS  DC 


The  machine  equations  used  in  the  computer  are 
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va  =  Vc 
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vp  =  2.09 

6X104  l*tgs 

[§] 
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The  analog  computer  diagram  of  equations  (APPII-15) 
through  (APPII-21)  is  presented  in  Figure  19  of  this  appendix. 
The  helicopter  voltage  response  (Vjj/2000)  was  obtained  to 
step  inputs  of  the  natural  charging  current  (In/5  x  10*5). 

As  expected*  the  system  stability  increased  with  a 
decrease  of  tg.  A  value  of  tg  ^  0.1  second  provides  good 
system  stability. 
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